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ABSTRACT 

The  gas  phase  hydration  of  the  alkali  ions,  M+ 

(M  =  Li,  Na,  Rb ,  Cs)  was  studied  by  means  of  a  specially 
designed  mass  spectrometer  which  was  capable  of  operating 
at  ion  source  pressures  in  the  torr  range.  The  mass 
spectrometer  was  used  to  analyze  the  equilibrium  concent¬ 
rations  of  the  gas  phase  hydrated  ions  arising  from  the 
thermal  equilibrium  of  alkali  ions  and  water  vapour  in 
the  ion  source.  The  analyses  proved  the  occurrence  of 
the  general  reaction  (1) 

(1)  M+(H20!n-l  +  H2°  *  M+(H2°>n  n  =  1-6 

where  notation  n-l,n  is  used  to  denote  specific  reactions. 

Equilibrium  constants,  Kn_-^  n ,  for  the  general  reaction 

(1)  were  determined  over  a  range  of  water  pressures. 

Typically  the  pressure  was  varied  from  0.2  to  3  torr  to 

demonstrate  that  n  is  pressure  independent.  By  the 

determination  of  the  equilibrium  constants  K  ,  at  various 

^  n-l,n 

O  O 

temperatures  the  thermodynamic  values  AH^_^  ,  AGn_^  n  and 

O 

AS  n  were  obtained. 

n-1,  n 

A  theoretical  calculation  of  the  potential  energies 
+ 

of  the  clusters  M  (H^O) r  is  presented.  This  calculation 
was  based  on  the  various  electrostatic  forces  between  the 
central  ion  and  the  ligand  molecules  as  well  as  those 

O 

among  the  ligand  molecules.  The  potential  energies  AEn_-^  n 

o 

were  compared  with  AH^_^  ^  obtained  experimentally.  The 


Ill 


o  o 

comparison  of  A.Hn_^  ^  with  AE^_^  ^  revealed  some  important 
features  concerning  the  nature  of  hydration  of  alkali  ions 
in  the  gas  phase,  i.e.  the  type  of  bonding  and  the  struct- 
ure  of  M  CH2°)n  clusters  where  n  =  1  to  6.  The  entropy 
values,  ASq  ^  for  reaction  (1)  calculated  from  statistical 
mechanics  were  found  to  be  in  excellent  agreement  with  the 
experimental  values. 

From  the  results  of  this  work  and  the  study  of  the 
gas  phase  hydration  of  halide  negative  ions,  carried  out 
by  other  workers  in  this  laboratory,  the  hydration  energies 
of  the  positive  and  negative  ions  were  compared.  Our 
studies  indicate  that  initial  hydration  steps  are  more 
favorable  for  positive  ions.  However  as  more  water  mole¬ 
cules  are  added  to  the  clusters,  negative  ion  hydration 
becomes  more  favorable. 

Another  important  feature  of  the  present  study  was 
the  comparison  of  the  gas  phase  results  for  the  hydration 
of  alkali  and  halide  ions  with  the  individual  ion 
hydration  energies  in  solution.  The  trend  in  our  data 
supports  the  individual  ion  hydration  energies  obtained  by 


Randles . 


IV 


ACKNOWLEDGEMENTS 


Sincere  appreciation  and  gratitude  to  Dr.  P.  Kebarle, 
whose  patience,  originality,  and  thoroughness  have 
been  the  guidance  and  encouragement  throughout  this 
research . 

The  author  is  deeply  indebted  to  Dr.  S.  K.  Searles 
who  schooled  the  author  in  the  technique  of  mass 
spectrometry . 

The  author  wishes  to  thank  Dr.  R.  Yamdagni  for  pro¬ 
gramming  the  theoretical  calculations  and  other 
members  of  the  mass  spectrometry  group  for  their 
valuable  suggestions  during  the  course  of  this  study. 

Special  thanks  go  to  Mrs.  Mary  Waters .  Her  care  in 
typing  the  manuscript  was  invaluable. 

The  financial  assistance  provided  by  the  University 
of  Alberta  is  gratefully  acknowledged. 


, 


V 


TABLE  OF  CONTENTS 

Page 

Abstract .  ii 

Acknowledgements  .  iv 

List  of  Figures .  ix 

List  of  Tables .  xi 

1.  INTRODUCTION 

1.1  The  present  study .  ^ 

1.2  The  nature  of  solvation .  1 

A.  Solvation  of  the  ions  in  the  solution.  i 

B.  Solvation  of  the  ions  in  the  gas  phase  4 

1.3  Principles  of  the  mass  spectrometric 

method  for  the  present  gas  phase  sol¬ 
vation  studies . 11 

1.4  Thermodynamic  relations  relevant  to 

gas  phase  solvation  reactions  ....  13 

1.5  Ion-solvent  molecule  interaction 

studies  in  the  present  laboratory  .  .  14 

+ 

1.6  Purposes  of  the  M  n  study  .... 

1.7  Significance  of  gas  phase  solvation 

reactions . . 19 

2 .  EXPERIMENTAL 

2 . 1  General . 24 

2.2  Instrumental  requirements  to  observe 

gas  phase  hydration  equilibria  and  the 
principles  of  method  .  25 

2.3  The  reaction  chamber . 30 


VI 


Page 

2.4  The  reaction  chamber  pressure  measure¬ 
ment  . 33 

2.5  The  reaction  chamber  temperature  con¬ 
trol  and  measurement . 33 

2.6  Creation  of  alkali  ions  and  manufact¬ 
ure  of  a  thermionic  source . 35 

2.7  Production  and  transport  of  the  hyd¬ 

rated  ions  out  of  the  reaction 

chamber . 37 

2.8  Ion  acceleration,  mass  analysis  and 

ion  current  detection  .  41 

2.9  Relation  of  the  ion  current  to  the 

ion  concentration . 42 

A.  Discrimination  in  transport  through 

the  reaction  chamber  .  42 

B.  Discrimination  due  to  stripping  out¬ 
side  of  the  reaction  chamber . 43 

C.  Discrimination  in  mass  analyzer  tube  .  44 

D.  Discrimination  in  the  ion  current 

detection . 45 

3.  RESULTS  OF  HYDRATION  OF  POSITIVE  ALKALI 

IONS  M+ . 47 

3.1  Introduction  and  thermodynamic 

relations . 47 

3.2  Presentation  of  the  results  and 

general  discussion  .  48 

3.3  Equilibrium  distribution  of  M^CfL^O)  .  83 


' 

i  1 ' 

. 

Vll 


Page 


4.  ELECTROSTATIC  CALCULATIONS  OF  THE  POTENTIAL 

ENERGIES  OF  THE  CLUSTERS  M+ (Ho0)  AND  COM- 

2  n 


PARISON  WITH  EXPERIMENTAL  RESULTS  .  86 

4.1  Introduction . 85 

4.2  Interaction  energies  .  85 

A.  Ion-permanent  dipole  forces  .  86 

B.  Ion-induced  dipole  forces  .  86 

C.  Dispersion  forces  .  87 

D.  Dipole-dipole  repulsive  forces  .  87 


E.  Ion-ligand  molecule  repulsive  forces  .  .  87 

4.3  The  present  model  for  the  electrostatic 

calculations . 88 

A.  Ion-permanent  dipole  interaction  (EDIP)  .  88 

B.  Ion-induced  dipole  interaction  (EPOL)  .  .  91 

C.  Ion-ligand  and  ligand-ligand  interaction 

energy  due  to  dispersion  forces  (EDIS) .  .  92 

D.  Dipole-dipole  repulsion  energy  (RDIP)  .  .  93 

E.  Ion-ligand  molecule  repulsive  energies 

(REL) . 96 

4.4  Evaluation  of  the  numerical  value  for  the 
constant  A  in  the  electronic  repulsion 

term  REL  =  AR“12 . 9  7 

4.5  Discussion  of  electrostatic  calculations . 10 4 

4.6  Comparison  of  electrostatic  potential 

energies  with  previous  calculations  .  .  .110 

4.7  The  comparison  of  calculated  E°  ,  with 

o  +  n  r ,  n 

experimental  AH  ,  for  M  (H„0)  ,  + 

n  _L  j  n  ^  n  -L 

H20  X  M  (H20)n 


112 


' 


Vlll 


Page 

4.8  Calculation  of  the  entropy  of  the  M+(H20)  .  120 

5.  APPLICATION  OF  EXPERIMENTAL  RESULTS  TO  PROCESSES 

IN  AQUEOUS  SOLUTION  .  130 

5.1  Comparison  between  the  hydration  energies 

of  the  alkali  and  halide  ions . 130 

5.1a  in  solution . 130 

5.1b  in  the  gas  phase . 131 

5.2  Ion-ligand  bond  energies  and  hydration 

numbers . 140 

6.  OTHER  APPLICATIONS  OF  HYDRATION  MEASUREMENTS  IN 

THE  GAS  PHASE . 146 

6.1  Clustering  equilibria  and  the  catalysis  of 

water  condensation  by  ionic  charges  ....  146 

6.2  Proton  affinities  of  the  alkali  hydroxides 

in  the  gas  phase . 153 

7.  SUGGESTION  FOR  FURTHER  RESEARCH  .  159 

7.1  Solvation  of  alkali  and  halide  ions  by 

polar  aprotic  solvents . 159 


BIBLIOGRAPHY 


161 


' 

Figure 


LIST 


0  F 


FIGURES 


2.1 

2.2 

2.3 


3.1 

to 

3.21 


3 .22 

to 

3.25 

3.26 


3.27 


3.28 


3.29 


3.30 


4.1 

4.2 


Page 

Cross-sectional  view  of  mass  spectrometer.  26 

Reaction  chamber  and  acceleration  plates 

for  hydration  reactions  of  alkali  ions.  28 


Cross-section,  top  and  bottom  view  of  reaction 
chamber  32 


Log  K  ,  for  the  gas  phase  hydration  of 
n  l  r  n 

alkali  ions,  M  (M  =  Li,  Na,  Rb,  Cs),  at 
various  temperatures,  versus  pressure  of 
h20 

Van ' t  Hoff  type  plots  of  the  equilibrium 
constants,  Kn_^  n  for  gas  phase  hydration 

of  alkali  ions,  M+  (M  =  Li,  Na,  Rb,  Cs) 

Enthalpies  for  clustering  reactions 

M+(H20)n-l  +  H2°  *  M+(H2°)n 

Free  energies  of  clustering  reactions 

M+(Ho0)  .  +  H„0  •*  M+(Ho0) 

2  n-1  2  2  n 

Plot  of  AH  (Cs+)  -  AH  (M+)  versus  n 

o ,  n  o  ,  n 

showing  that  the  above  differences  approach 
differences  of  total  single  ion  heats  of 
hydration  obtained  by  Randles  (R)  and 
Latimer  (L) 

Plot  of  enthalpies  for  clustering  reactions 

M+(H„0)  +  H„0  -*■  M+(Ho0)  versus  AH  total 

v  2  n-1  2  2  n 

single  ion  heats  of  hydration  obtained  by 

Randles . 

Plot  of  equilibrium  distribution  of 
Na+(H20)n  at  T  =  300°K  as  a  function  of 
water  pressure 

M+(H20)  model  used  in  the  theoretical 
calculations . 

Plot  of  A2  on  logarithmic  scale  versus 
the  Van  der  Waals  radius  of  the  inert  gas 
atoms . 


49-69 


70-73 

79 

80 


81 


82 


84 

90 


100 


' 


X 


Figure 


Page 


4.3 

4.4 

4.5 

4.6 

4.7 

5.1 

5.2 

а, b 

б. 1 


6.2 

6.3 


Comparison  between  experimental  enthalpies 
AH  and  potential  energies  AE  obtained 

on  the  basis  of  electrostatic  type  calcu¬ 
lations  for  sodium  hydrates  Na+(H20)n 

Comparison  between  experimental  enthalpies 

AH  and  potential  energies  AE  for 
o,n  +  o,n 

cesium  hydrates  Cs 

Plot  of  potential  energies  evolved  for  the 
reactions:  M+  +  H^O  ->  M  (H^O) 

O 

Comparison  of  experimental  AH^_^  n  and 

calculated  potential  energy  differences 

°  + 

AE  -»  for  alkali  ion  clusters  M  (H~0) 
n-1 ,  n  2  n 

Principal  axes  of  H^O  and  Li+(H20)  used  to 
calculate  the  moments  of  inertia. 

O 

Comparison  of  n  f°r  hydration  of 

alkali  positive  and  halide  negative  ions 

Comparison  of  AHq  ^  hydration  of  alkali 

and  halide  ions  with  total  enthalpies  of 
hydration  of  Latimer  (L)  and  Randles  (R) 

Comparison  of  free  energy  of  hydration 

O 

AG  with  free  energy  of  evaporation  of 

n-l,n 

1  mole  of  water  molecules  from  one  mole  of 
water  droplets  carrying  a  single  ionic 
charge . 

The  variation  of  r/  00  as  a  function  of 
droplet  radius 

Proton  affinities  of  MOH  versus  ionization 
potential  of  M 


108 

109 

116 

117 

123 

132 

137-138 

147 

151 


156 


' 


XI 


Table 


3.1 
to 
3.4 

4.1 


4.2 


4.3 


4.4 


4.5 


4.6 


4.7 


4.8 


LIST  OF  TABLES 


Thermodynamic  values  for  the  gas  phase 
hydration  of  the  alkali  ions 


Ionization  potentials  and  polarizabilities 
of  alkali  ions 


Coefficient  "A"  for  repulsion  term  AR 
between  two  atoms  of  the  noble  gases 

Potential  energies  of  clusters  from 

electrostatic  calculations  using  repulsion 
term  AR  ^ 


Individual  energy  contributions  to  potential 
energy  calculations  of  cluster  Na^t^O) 

+ 

Potential  energies  of  clusters  M  (H2^n  ^rcm 
electrostatic  calculations  using  repulsion 
term  AR  ^ 

Moments  of  inertia  and  M — equilibrium 
distances 


Force  constants,  vibrational  frequencies  and 
entropy  contributions  due  to  stretching 
vibrations  of  M+ - OH^ 


Potential  energy  barriers  for  rotation  (AP.E.) 
vibrational  frequencies  and  entropy  contri¬ 
butions  due  to  two  torsional  vibrations  of 


Page 

75-78 

94 

99 

102 

106 

115 

122 

125 


127 


. 


Xll 


Table  Page 


4.9  Translational,  rotational  and  vibrational 

contributions  to  experimentally  determined 

°  +  + 

AS298°K  f°r  reacti°n:  M  +  HO  M  (OH  )  128 

5.1  Ion-water  molecule  bond  energies  in  the 

liquid  phase  143 

5.2  Hydration  numbers  "h "  for  alkali  ions  144 

o  o 

6.1  Free  energy  changes  AG  ^  ,  AG^  and 

AG  as  a  function  of  radius  of  charged 

evap  ^ 

droplet  150 

6.2  Proton  affinities  of  alkali  hydroxides  MOH 
and  thermochemical  data  used  for  their 
evaluation 


155 


1. 


1.  INTRODUCTION 
1 .1  The  Present  Study 

The  present  work  deals  with  the  hydration  energies 
of  the  alkali  ions  in  the  reactions  (1.1  to  1.4)  occurring 
in  the  gas  phase  where  M  stands  for  Li,  Na,  Rb  and  Cs . 


rH 

• 

i — 1 

M+ 

+ 

h2o 

£ 

M+ (H20) 

i — 1 

o 

1.2 

m+(h2o) 

+ 

h2o 

-> 

m+(h2o)2 

(1,2) 

1.3 

m+(h20)2 

+ 

o 

<N 

K 

-> 

<- 

m+(h20)3 

(2,3) 

• 

I — 1 

M+(H2°)n_i 

+ 

h20 

-y 

M+ (H„0) 

2  n 

(n-1 , n) 

The  present  work  is  a  continuation  of  an  earlier  investi¬ 
gation  which  dealt  with  the  hydration  of  the  potassium  ion 
(1) .  The  equilibrium  constants  for  attachment  of  water 
molecules  to  the  alkali  ions  were  determined  as  a  function 
of  temperature.  The  enthalpy,  free  energy  and  entropy 
changes  of  the  hydration  reactions  were  then  obtained  by 
Van't  Hoff  plots.  The  experimentally  determined  enthalpy 
and  entropy  changes  are  compared  with  theoretically  calcu¬ 
lated  values.  The  enthalpies  of  hydration  (n-l,n)  are  also 
compared  with  the  total  enthalpies  of  hydration  of  ions  in 
aqueous  solution. 

1 . 2  The  Nature  of  Solvation 
A.  Solvation  of  ions  in  the  solution 

Hydration  in  solution  is  the  result  of  interactions 
between  the  solute  (the  ions)  and  the  solvent  (water).  Not 
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only  must  the  interactions  between  the  ions  and  the  water 
molecules  in  their  immediate  vicinity  be  considered,  but 
also  the  interactions  with  all  the  rest  of  the  water.  Sol¬ 
vation  (hydration)  in  solution  may  be  understood  as  the 
sum  of  all  changes  caused  by  the  appearance  of  the  ions  of 
the  electrolyte  in  solution  (2) . 

The  knowledge  of  the  existence  of  hydrated  ions  in 
solution  dates  back  to  the  beginning  of  this  century  (3) . 
Despite  the  general  acceptance  that  ions  in  solution  are 
hydrated,  serious  disagreement  between  authors  arises  as 
soon  as  quantitative  treatments  of  hydration  are  attempted. 
This  difficulty  occurs,  because  the  water  molecules  in  the 
hydration  shell  are  in  continuous  exchange  with  the  water 
in  the  bulk  of  the  solution  and  it  is  only  on  a  time  average 
that  a  hydration  shell  may  be  said  to  exist  around  most 
ions  (4) .  In  order  to  account  for  the  large  discrepancy 
between  the  different  methods  of  measuring  hydration  in 
solutions,  Bockris  (5)  has  suggested  that  hydration  be 
divided  into  two  main  types,  which  he  called  "primary 
hydration"  and  "secondary  hydration".  "Primary"  or 
"close"  hydration  refers  to  the  water  molecules,  which  are 
immobilized  by  the  ion  and  move  as  one  entity  with  it, 
while  "secondary"  or  "distant"  hydration  refers  to  other 
ion-water  interactions  beyond  the  primary  hydration  shell. 

"Distant"  hydration  does  not  depend  on  specific 
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molecular  interactions  between  the  ion  and  the  solvent 
molecules  and  it  can  usually  be  dealt  with  satisfactorily 
by  classical  electrostatic  theories,  such  as  the  Born 
equation  (6),  which  assumes  a  continuous  solvent.  Theories 
based  on  continuous  solvent  models  must  cope  with  effects 
such  as  dielectric  saturation,  the  true  radii  of  ions 
in  solution  and  the  lack  of  spherical  symmetry  near  the 
surface  of  the  ion  arising  from  the  finite  size  of  the 
water  molecules  (6) . 

On  the  other  hand,  "close"  hydration  does  depend  on 
specific  molecular  interactions  between  the  ion  and  the  sol¬ 
vent  molecules.  Molecular  interaction  models  must  take 
into  account  the  coordination  number  and  the  orientation  of 
the  water  molecules,  ion-dipole,  dipole-dipole  inter¬ 
actions  etc  (6) . 

For  the  over-all  effect  of  hydration  Bockris  used 
the  expression  "total  hydration". 

The  separation  of  the  "total  hydration"  effect  into 
the  interactions  between  ions  and  water  molecules  in  their 
immediate  vicinity,  and  the  interactions  with  the  rest  of 
the  water  molecules  represents  a  simplification  (6) .  Since, 
however,  in  liquids  increasing  distance  between  two  given 
particles  is  attended  by  a  continuously,  but  gradually, 
decreasing  regularity  some  simplification  must  be  made  to 
define  a  boundary  between  the  water  molecules  of  "close" 
and  "distant"  hydration. 
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It  will  be  shown  in  the  next  section  that  studies 
of  the  hydration  of  ions  in  the  gas  phase  allow  determin¬ 


ations  of  the  thermochemistry  and  reveal  a  relative 
stability  of  "close"  hydrated  ions  in  the  absence  of  the 
bulk  solvent. 

B .  Solvation  of  ions  in  the  gas  phase 

The  interactions  of  ions  with  solvating  molecules 
in  solution  have  been  studied  for  many  years.  The  cor¬ 
responding  studies  of  interaction  of  gaseous  ions  with 
"solvent"  vapor  are  relatively  recent  (7). 

Ions  composed  of  clusters  or  a  group  of  molecules 
of  the  type  M+(H20)n  where  n  is  a  small  number,  must 
exist  in  the  gas  phase  due  to  the  forces  of  attraction 
between  ions  and  neutral  molecules.  An  ion  of  charge  e 

has  associated  with  it  an  electrical  field  of  strength 
2 

F  =  e/r  at  a  distance  r.  This  field  is  of  the  order  of 

o  _8  4  —6 

10  volts/cm  at  r  =  10  cm  and  10  volts/cm  at  10  cm. 

In  fields  of  this  magnitude  all  molecules  experience  a 
displacement  of  their  electron  cloud  relative  to  the 
nuclei  so  as  to  produce  an  induced  electrical  dipole 
moment  of  pi  =  el  =  otF,  where  1  represents  the  effective 
separation  of  the  positive  and  negative  center  and  a 
the  polarizability  of  the  molecule.  In  addition  a  mole¬ 
cule  like  H20  will  also  have  an  effective  part  of  its 
permanent  dipole  moment  yp,  oriented  by  the  electric  field 
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of  the  ion.  Thus,  between  the  ion  and  molecule,  there  is  a 
force  of  attraction  due  to  ion-permanent  dipole  and  ion- 
induced  dipole  forces . 

The  experimental  indications  for  the  existence  of 
clustered  ions  were  obtained  by  the  earliest  workers  in 
the  field  of  gaseous  electronics.  At  the  beginning  of 
this  century  Rutherford  and  Thomson  (8)  discovered  that 
the  mobilities  of  ions  in  gases  were  lower  than  values  pre¬ 
dicted  by  the  kinetic  theory.  They  attributed  this 
deviation  to  the  formation  of  ion-clusters.  The  formation 
of  ion-clusters  was  explained  on  the  basis  of  the  attrac¬ 
tion  forces  discussed  above. 

The  rough  criterion  for  the  stability  of  these  ion- 

clusters  was  set  by  the  earliest  workers  (10)  as  follows: 

PE 

——  =  1.  PE  represents  the  potential  energy  of  a  molecule 

in  the  cluster  and  KE  stands  for  the  kinetic  energy  of 

average  ion  molecule  impact.  This  comes  from  the  Boltzman 
AE 

factor  e  KT  ,  where  AE  =  -PE  and  RT  =  KE .  In  essence  the 
above  expression  states  that  if  PE  =  KE  then  there  are 
roughly  equal  numbers  of  clustered  and  free  ions. 

Already  in  the  early  1920's  it  was  recognized  that 
the  above  criterion  was  faulty  (9) .  Thermodynamic  equil¬ 
ibrium  between  an  ion  cluster  and  its  dissociation  products 
essentially  depends  not  only  upon  the  Boltzman  factor  but 
also  upon  entropy  factors  which  weight  it  heavily  in 
favour  of  dissociation.  This  means  that  the  a  priori 
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probability  of  a  state  must  be  included  as  a  factor  in  the 
ratio  of  PE/KE  in  order  that  the  condition  for  stability 

be  fully  represented  (9) .  The  ion  clusters  are  a  "condensed" 
state  with  relatively  low  entropy,  while  the  dissociated 
ions  and  free  molecules  constitute  a  "dispersed"  state 
with  a  higher  entropy,  or  statistical  weight.  Unless  much 
is  known  of  the  specific  nature  of  a  given  possible  cluster 
and  its  constituents,  the  number  of  molecules  etc.,  it  is 
impossible  to  evaluate  the  relative  statistical  weight 
factor  g  for  dissociated  and  condensed  states.  It  is 
accordingly  necessary  to  set  the  criterion  for  cluster 
stability  as  PE/gKE  _>  1  (9)  .  In  this  relation,  the 
relative  statistical  weight  g  of  the  dispersed  state  may 
be  greater  than  one  by  a  considerable  margin. 

A  second  theory,  to  explain  the  discrepancy  between 
the  theoretical  and  experimental  values  of  ionic  mobility 
was  put  forward  by  Sutherland  and  Wellisch  (9) .  This 
theory,  which  later  on  became  known  as  "small  ion  theory", 
assumes  that  the  increased  size  of  the  clustered  ions  is 
not  the  major  factor  in  retarding  the  motion  through  the 
gas,  but  that  the  motion  is  retarded  by  momentum  transfer 
as  a  result  of  the  electrical  forces  between  the  ions  and 
neutral  molecules. 

Munson,  Tyndall  and  Hoselitz  (10)  studied  ionic 
mobility  in  which  the  gas  phase  transport  of  ions  was 
investigated  as  a  function  of  F/P;  F  is  the  electric  field 
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strength  and  P  is  the  gas  pressure.  For  example,  in  the 

study  of  the  mobility  of  alkali  metal  ions  in  "moist" 

rare  gases,  the  ion-clusters,  M+(H„0)  ,  n  <  6  were  used 

n 

to  explain  the  observed  ionic  mobilities.  Munson  et  al . 
observed  that  the  mobility  of  the  clustered  ion  is  not 
constant,  but  shows  a  variation  with  F/P.  At  high  F/P 
only  the  unclustered  ions  appear  and  at  low  F/P  only 
the  "fully  clustered"  ions  appear.  To  explain  the  anomaly 
of  the  appearance  of  only  two  mobilities  and  yet  the 
variation  of  the  mobility  with  F/P  the  following  hypothesis 
was  made  (10) .  A  cluster  cannot  form  readily.  It  is 
assumed  that  until  the  first  molecule  is  attached  to  the 
ion,  there  will  be  no  attachment  for  the  other  molecules. 
When  the  first  molecule  attaches  to  the  alkali  ion,  it 
must  free  itself  from  the  energy  of  binding.  Therefore  a 
third  body  impact  is  required  to  release  this  energy,  other¬ 
wise  the  ion-cluster  dissociates.  Once  the  first  molecule 
has  attached  itself  the  other  molecules  can  also  do  so; 
the  excess  binding  energy  can  be  taken  up  by  the  cluster 
as  internal  vibrational  energy,  which  can  be  transferred 
to  surrounding  molecules  in  later  collisions.  From  the 
initial  addition  of  the  first  H20  molecule  in  a  three- 
body  impact,  where  a  noble  gas  atom  serves  as  third  body, 

the  other  f^O  molecules  in  the  cluster  can  attach  as 
rapidly  as  encountered  up  to  the  equilibrium  number.  The 
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average  number  of  these  molecules  in  the  resulting  labile 
cluster  will  vary  as  the  energy  of  the  ion  varies.  Below 
the  highest  F/P,  possibly  only  the  original  ion  will  remain 
with  one  labile-attached  molecule.  These  will  be  attached 
and  detached  in  rapid  sequence,  giving  a  continuously  vary¬ 
ing  mean  mobility  of  the  clustered  ions  for  each  F/P  value. 
Therefore,  the  change  of  the  mobility  of  clustered  ions 
with  F/P  according  to  Munson  et  al .  (10)  is  due  to  changing 

the  average  size  of  the  cluster  formed. 

The  average  size  of  the  ion-cluster  varies  not  only 
with  the  polarity  of  the  surrounding  molecules  but  also 
with  the  size  of  the  ion  and  the  temperature.  This  raises 
the  question  as  to  whether  only  strongly  polar  molecules 
can  form  ion-clusters.  Some  light  was  thrown  on  this 
matter  by  the  work  of  Munson  and  Hoselitz  (11) ,  who  studied 
the  mobilities  of  alkali  ions  at  low  fields  F,  in  various 
nonpolar  gases.  The  presence  of  clustering  was  observed 
at  low  temperatures  in  some  gases  and  even  at  room 
temperatures  with  Li  ions  in  Xe  gas . 

Although  the  existence  and  the  importance  of  ion 
clusters  in  mobility  measurements  has  been  recognized  for 
many  years,  the  systematic  study  of  ion-molecule  clusters 
in  the  gas  phase  did  not  start  until  development  of  "high 
pressure"  mass  spectrometry.  In  an  ordinary  mass  spect¬ 
rometer  used  for  analysis  (isotope  analysis,  qualitative 
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and  quantitative  analysis  of  vapours)  the  pressure  in  the 
ion  source  and  mass  analysis  region  is  kept  low  so  that 
the  ions  can  leave  the  ion  source  and  be  mass  analyzed 
without  the  occurrence  of  collisions  with  gas  molecules. 

If  the  pressure  of  the  ion  source  is  increased,  collisions 
of  the  ions  with  gas  molecules  begin  to  occur  and  these 
can  lead  to  ion  molecule  reactions  which  change  the 
nature  of  the  primary  ions.  The  first  ion-molecule 
reaction  was  observed  by  J.  J.  Thomson  (12)  with  the 
first  mass  spectrograph.  However  it  was  only  recognized 
during  the  early  1950's  that  ion-molecule  reactions 
might  play  a  large  role  in  gas  phase  chemistry.  With 
this  in  mind  V.  Talroze  (13),  D.  P.  Stevenson  (14), 

W.  H ami 11  (15)  and  F.  H.  Field,  J.  L.  Franklin  and  F.  W. 
Lampe  (16)  started  systematic  studies. 

For  the  purposes  of  the  present  discussion  ion- 
molecule  reactions  can  be  divided  into  two  types:  true 
chemical  reactions  in  which  chemical  bonds  are  formed 
and  broken  and  clustering  reactions  where  association 
occurs  mainly  because  of  the  operation  of  electrostatic 
forces  i.e.  ion-dipole  and  ion-induced  dipole  attractions. 
An  example  of  a  chemical  process  is  given  by  the  associ¬ 
ation  of  an  ethylene  ion  with  an  ethylene  molecule  to 
form  a  butene  ion  (reaction  1.5),  while  reactions  1.6 
and  1.7  represent  typical  clustering  reactions 
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1.5 

C2H4+ 

+ 

C2H4 

ll 

n 

rc 

+ 

8 

1.6 

Li+ 

+ 

nH20 

=  Li+ 

(H20) 

1.7 

Li  + 

+ 

nXe 

=  Li  + 

(Xe)  n 

It  is  interesting  to  note  that  almost  all  ion-molecule 
investigations  have  dealt  with  reactions  of  the  chemical 
type. 


The  clustered  ions  Ho0+(Ho0)  have  been  detected 

3  2  n 

mass  spectrometrically  by  Beckey  (17)  who  observed 
clusters  H30+(H20)n  up  to  n  =  10.  Knewstubb  and  Sugden 
(18)  in  a  mass  spectrometr ic  study  of  ionization  in  flames 

-f-  4" 

observed  the  ions  H_,0  (H~0)  and  NH.  (NH~)  up  to  n  =  4 . 

3  2  n  4  3  n  r 

Sieck,  Abramson  and  Futrell  (19)  have  studied  the 

+ 

mechanisms  leading  to  the  formation  of  H  (CH^OH)  and 

H+ (CnHcOH)  - .  The  ions  H+(H„0)  ,  H+ (CHo0H)  ,  H+(CHo0CHo)  , 

zoo  z  n  on  o  o  n 

H+ (CH_NH0 )  ,  H+[(CH_)„NH]  and  H+[(CH_)^N]  were  observed 

v  3  27  n'  3  2  n  3  3  n 

by  Munson  (20)  in  the  high  pressure  mass  spectra  of 
CH3OH,  CH3OCH3,  CH3NH2,  (CH3)2NH  and  (CH3)3N  respectively, 
where  n  =  1-2  for  each  species.  Yang  and  Conway  (21)  have 
studied  the  system  °2n+2  up  to  n  =  4 .  The  clustered  ion 
K+(H20)  was  detected  by  Chupka  (22)  in  a  mass  spectro- 
metric  study  of  ions  emerging  from  a  heated  Knudsen  cell 
containing  KC1  salt.  The  ion-clusters  M+(H20),  where 
M  =  Li,  Na  and  K  were  observed  in  flames  by  Hayhurst  and 
Sugden  (23) .  They  examined  the  extent  to  which  alkali 
ions  are  hydrated  under  flame  conditions  (premixed  H2,C>2, 
N2)  at  atmospheric  pressure.  Aqueous  solutions  of  alkali 
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metal  chlorides  were  sprayed  into  the  gas  mixture  to  give 

a  concentration  of  alkali  chloride  in  the  flame  of  1  part 
7 

in  10  .  The  principal  positive  ions  detected  by  flame 

mass  spectrometry  were  (H20) .  From  the  constancy  of 

the  ratio  M+/M  (H^O)  with  distance,  over  a  restricted 

range  of  distances,  near  the  reaction  zone  in  which  zone 

the  temperature  was  nearly  constant,  it  was  assumed  that 

the  equilibrium  for  the  reaction  M+  +  H^O  ■+  M+(H20)  was 

established.  On  this  assumption  equilibrium  constants, 

K  were  derived  for  two  temperatures.  From  the  vari- 
eq  L 

ation  of  equilibrium  constants,  Kn  ,  with  temperature 

0/1 

AHq  ^  values  of  54,  27  and  40  kcal/mole  were  obtained 

-f  ,  , 

for  M  =  Li ,  Na,  K  respectively.  The  error  in  the 

above  values  for  AH^  ^  is  apparent,  since  the 
+  + 

AHq  ^(K  )  >>  AHq  ^ (Na  ) .  This  is  opposite  to  what  one 

would  expect,  since  Na+  is  a  smaller  ion  than  K+  and 
therefore  Na+  -  H^O  interaction  energy  should  be  bigger 
than  K+  -  H^O  interaction.  This  indicates  that  the 
author's  assumptions  about  equilibrium  conditions  and 
the  constancy  of  the  temperature  in  the  reaction  zone 
are  somewhat  illusory. 

1 . 3  Principles  of  the  Mass  Spectrometric  Method  for 

the  Present  Gas  Phase  Solvation  Studies 
Mass  spectrometric  observations  of  the  clustered 
ions  discussed  in  the  previous  section  were  mainly  in- 
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cidental  and  were  not  pursued  further.  The  most  exten¬ 
sive  studies  of  the  ion-cluster  or  ion  solvation  reactions 
in  the  gas  phase  occurring  under  equilibrium  conditions 
were  carried  out  in  this  laboratory. 

In  order  to  measure  the  gas  phase  solvation  of 
ions,  the  instrument  required  must  produce  ions,  allow  them 
to  react  and  equilibrate  at  high  pressures  and  finally 
analyze  the  reactant  and  product  ions.  Because  of  its 
ability  to  satisfy  all  three  requirements  a  modified  mass 
spectrometer  has  been  used  in  this  study. 

The  alkali  ions  were  produced  by  thermionic  emission 
from  a  heated  filament  coated  with  suitable  salt  containing 
the  desired  ion.  Since  thermionic  emission  can  be 
obtained  only  at  elevated  temperatures,  the  thermionic 
source  was  separated  from  the  reaction  chamber  whose 
temperature  was  controlled  independently.  In  order  to 
attract  the  ions  from  the  thermionic  source  to  the 
reaction  chamber,  a  potential  difference  was  applied  between 
thermionic  source  and  reaction  chamber.  The  ions  entering 
the  field  free  reaction  chamber  gradually  lose  their  excess 
kinetic  energy,  through  collisions  with  neutral  water 
molecules.  Thermalized  ions  diffusing  to  the  exit  slit 
and  escaping  into  the  vacuum  region  are  captured  by 
electric  fields,  accelerated,  and  magnetically  mass  analyzed. 
In  the  present  measurements  the  pressures  used  in  the 
reaction  chamber  were  in  the  range  of  0.2  -  3  torr .  Since 
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the  ion  acceleration,  mass  analysis  and  collector  regions 

-4 

should  be  maintained  at  less  than  10  torr  for  the  mass 
spectrometer  to  function  properly,  the  present  instrument 
used  a  high  capacity  pumping  system  and  reaction  chamber 
with  a  narrow  exit  slit  (to  reduce  the  gas  flow  out  of  the 
reaction  chamber) . 


1 . 4  Thermodynamic  Relations  Relevant  to  Gas  Phase  Solvation 

Reactions 

Alkali  ions  which  are  ejected  by  a  heated  filament 
drift  into  the  reaction  chamber  under  the  influence  of  an 
electric  field.  In  the  field  free  reaction  chamber  the 
alkali  ions  lose  their  excess  translational  energy  and 
react  with  water  molecules  to  form  hydrated  ions.  Assuming 
that  equilibrium  is  achieved,  the  equilibrium  constant 
Kn-1  n  ^or  t^ie  hydration  reaction  n-l,n  (reaction  1.4) 
is  given  by  the  following  relation: 


K 


n-1 ,  n 


PM+ (Ho0) 

2  n 

p  +  P 

M  (Ho0)  ,  Ho0 

2  n-1  2 


is  the  pressure  of  the  reacting  ion  and 
that  of  the  product.  is  the  partial  press 

of  water  in  the  reaction  chamber.  The  pressure  ratio  of 


where  PM+(H  0) 

2  n-1 


PM+ (Ho0) 

2  n 


the  ions  was  assumed  to  be  equal  to  the  ratio  of  the 


intensities  of  the  respective  ions,  recorded  by  the  mass 
spectrometer.  The  validity  of  this  assumption  is  further 


, 
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discussed  in  Chapter  2. 

The  Gibb's  free  energy  of  the  reactions  1.1  to  1.4  can 
be  calculated  by  the  following  equation: 


(1.1) 


AG 


n-1 ,  n 


-RT  log 


K 


n-1 ,  n 


where  R  is  the  universal  gas  constant,  T  is  the  absolute 


temperature  and  K  is  the  equilibrium  constant. 

r  n-l,n  ^ 

Since  the  values  for  K  n  in  the  present  work  were 

n-l,n  ^ 

calculated  for  a  standard  state  of  1  torr,  the  values  for 


K  *i  were  multiplied  by  760  before  substitution  into 

n-l,n  ^  J 

O 

the  above  equation  in  order  to  obtain  AG  .  in  a 
^  n-l,n 

standard  state  of  1  atm. 


From  the  measurements  of  the  equilibrium  constants 


at  various  temperatures,  it  was  possible  to  determine 


the  enthalpies  AH^_^  n  by  means  of  Van ' t  Hoff  plots.  The 

O 

entropy  values  AS  n  were  calculated  from  the  well  known 
^  n-l,n 


/ 1  9  \  AG  ,  =  AH  ,  -  TAS 

n-l,n  n-l,n  n-1, 


n 


1 . 5  Ion-solvent  Molecule  Interaction  Studies  Done  in  the 
Present  Laboratory. 

In  the  past  few  years  the  following  ion  solvation 
reactions  occurring  under  equilibrium  have  been  investi¬ 
gated  in  this  laboratory: 


H30+(H20)n 
The  H30+(H20)n 


system  has  been  extensively  studied  in 


' 


this  laboratory  (24)  and  thermodynamic  values  for  the 
reactions 
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have  been  determined  up  to  n=8.  More  recently  the 

H^0+ (H^O) ^  system  was  also  invesigated  by  F.  H.  Field  (25) 


and  his  results  are  in  excellent  agreement  with  Kebarle 
et  al  (24) .  A.  Good,  D.  Durden  and  P.  Kebarle  (26)  have 
made  a  kinetic  study  of  hydration  of  the  proton.  Employing  a 
quadrupole  mass  spectrometer  they  have  determined  the 
rate  constants  for  the  forward  reactions 


N 


+ 


H 


for  n=2  to  n=4 .  With  the  aid  of  the  equilibrium  constants 
determined  previously  (24)  in  this  laboratory,  they  were 
able  also  to  calculate  the  rate  constant  of  the  reverse 
reactions . 

In  addition  to  the  reactions  cited  above,  the  com¬ 
petitive  solvation  of  protons  by  H^O  and  CH^OH  was  also 
investigated  (27)  in  this  laboratory. 


laboratory  (28,29)  and  the  thermodynamic  values  for  the 


NH4+(NH3>n-l  +  NH3  t  NH4+(NH3>n 


reactions : 


' 


■ 
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have  been  determined  up  to  n=5.  In  addition  competitive 
solvation  of  NH^  +  by  NH^  and  H^O  was  also  investigated 
(28)  . 


The  gas  phase  hydration  of  halide  negative  ions  have 
been  studied  in  this  laboratory  by  M.  R.  Arshadi,  R. 
Yamdagni  and  P.  Kebarle  (30) .  By  determination  of  the 
equilibrium  constants  for  the  general  reaction: 


X  (H20) 


n-1 


+  H20  i  X  (H20)n 


at  various  temperatures  the  thermodynamic  values  of 

o  o  o 

AH  ,  ,  AG  ,  and  AS^  ,  were  determined.  In  addition 

n-l , n  n-l , n  n-l,n 

to  the  halide  negative  ions,  the  gas  phase  hydration  of 
C>2  ,  OH  and  NO^  were  also  studied  (31)  . 

1.6  Purposes  of  the  M+(H20)^  Study 

One  advantage  of  working  in  the  gas  phase  is  that 
the  hydration  of  positive  and  negative  ions  can  be 
studied  separately.  Although  positive  or  negative  ions 

have  a  separate  existence  in  the  gas  phase,  in  solution 
the  total  number  of  positive  charges  is  always  equal  to 
the  total  number  of  negative  charges  present.  This 
makes  it  very  difficult  to  assess  the  properties  of 
individual  ions  in  solution  since  one  is  always  dealing 
with  both  species.  The  question  of  the  possibility  of 


- 
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producing  only  positive  or  only  negative  ions  in  the  gas 
phase  i.e.,  the  separate  existence  of  the  ions  in  the  gas 
phase  is  due  to  very  low  concentrations  at  which  the  work  is 
possible.  Thus  it  can  be  calculated  that  the  ion  concen¬ 
trations  are  roughly  10^  ions/cc,  which  corresponds  to  a 
-14 

10  molar  concentration  in  solution.  At  such  low  con¬ 
centrations  also  separate  positive  ions  could  be  pro¬ 
duced  in  aqueous  solution,  for  example  by  shooting  in 
+ 

Na  ions.  However,  they  could  not  be  studied  in 
solution  at  these  low  concentrations.  The  advantage  of 
the  technique  thus  lies  in  the  ability  to  measure  extremely 
low  concentrations  of  ions  in  the  gas  phase  with  the  mass 
spectrometer . 

The  excellent  review  on  the  methods  of  assigning 
individual  ion  heats,  free  energies  and  entropies  of 
hydration  in  solution,  given  by  Rosseinsky  (31)  indicates 
that  most  methods  are  either  based  on  dividing  by  two 
the  thermodynamic  functions  of  selected  anion-cation 
pairs  or  on  using  some  extrapolation  procedure  involving 
the  ionic  radii.  Since  water  is  not  a  dimensionless  and 
structureless  liquid,  the  hydration  of  equal-size  cations 
and  anions  is  not  necessarily  the  same.  Some  authors  have 
attempted  to  account  for  the  specificity  in  ionic  hydration 
by  adding  different  increments  to  the  crystal  radii  of 
anions  and  cations  (32) ,  by  considering  the  different 
orientation  of  water  molecules  around  the  ions  (33) ,  or 
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by  taking  into  account  the  quadrupole  moment  of  the  water 
molecule  (34) . 

The  agreement  between  these  different  methods  is 
often  surprisingly  good,  especially  in  the  case  of  indiv¬ 
idual  heats  of  hydration.  However,  none  of  these  values 
should  be  considered  as  entirely  reliable  since  they  are 
all  based  on  the  assumption  that  the  true  ionic  radii 
are  known.  The  only  technique  that  has  been  developed 
for  the  evaluation  of  absolute  hydration  energies  of 
single  ions,  without  any  assumptions  about  the  dimensions 
of  the  ions,  is  based  on  the  measurement  of  the  volta 
potential.  The  volta  potential  is  measured  between  one 
electrolyte  solution  and  another  surface  separated  from 
the  solution  by  a  vacuum  or  an  inert  gas  (35) .  Randles 
(35)  has  applied  this  method  to  evaluate  individual 
hydration  energies  of  alkali  and  halide  ions. 

In  the  beginning  of  this  section  it  was  mentioned  that 
one  of  the  advantages  of  working  in  the  gas  phase  is  that 
the  hydration  of  positive  and  negative  ions  can  be  studied 
separately.  Thus  one  of  the  purposes  of  the  study  of 
the  hydration  of  the  positive  alkali  ions  was  to  compare 
the  results  with  the  hydration  energies  of  negative  halide 
ions.  The  gas  phase  hydration  of  halide  ions  was  carried 
out  in  this  laboratory  (30),  simultaneously  with  the 
present  study.  The  most  important  feature  of  the  present 
study  is  a  comparison  of  the  gas  phase  results  for  the 
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hydration  of  the  alkali  and  halide  ions  with  the  individual 
ion  hydration  energies  in  the  solution  evaluated  as  dis¬ 
cussed  above.  This  comparison  will  be  discussed  in 
chapter  5 . 

Another  important  feature  of  the  present  study  is  a 

O 

comparison  of  experimental  hydration  energies  AHn_^  ^  with 

O 

classical  electrostatic  energies  AE  ,  ,  based  on  ion- 

v  n-l,n' 

dipole  interactions.  Since  the  alkali  ions  are  spherical, 
one  might  expect  that  their  hydration  energies  should  be 
most  amenable  to  this  type  of  calculation. 

The  present  results  have  also  allowed  us  to  evaluate 
the  proton  affinities  of  the  alkali  hydroxides  in  the  gas 
phase,  discussed  in  chapter  6.  Proton  affinities  of 
alkali  hydroxides  have  not  been  previously  reported. 

The  applications  of  the  gas  phase  cluster  studies 
are  not  restricted  to  ionic  solutions.  There  are  a  number 
of  fields  in  which  ion  clusters  play  an  important  part. 
These  fields  are  considered  in  the  next  section. 

1 . 7  Significance  of  Gas  Phase  Solvation  Reactions 

Before  leaving  this  chapter,  it  is  appropriate  to 
mention  in  a  few  words  the  importance  of  ion  cluster 
studies  and  their  application  in  the  other  fields  of 
chemistry  and  physics.  The  role  of  clusters  in  ion 
mobility  and  flame  chemistry  has  already  been  mentioned. 
Lind  (36)  has  briefly  reviewed  the  role  of  ionic  clusters 
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in  radiation  chemistry.  According  to  Lind's  theory  an  ion 
can  attract  one  or  two  molecules  through  polarization 
and  hold  them  until  neutralized  by  an  oppositely  charged 
ion.  The  heat  of  neutralization  and  the  close  proximity 
of  ions  and  attached  molecules  would  promote  chemical 
reaction,  polymerization,  decomposition  etc.  Later  on 
this  theory  was  over-shadowed  by  the  free  radical  theory 
of  Eyring,  Hirschfelder  and  Taylor  (37) .  The  Eyring 
theory  postulates  that  the  reaction  of  single  ions  and  free 
radicals  and  molecules  is  more  probable  than  large  ion 
cluster  reactions.  The  product  of  the  neutralization 
of  ions  or  ion  clusters  and  also  the  products  of  dis¬ 
sociation  of  excited  molecules  may  be  free  radicals  or 
atoms  which  can  initiate  chain  reactions  and  thus  improve 
yield  (36 ) . 

Early  in  the  1950's,  however,  Lind's  theory  gained 
additional  support  when  it  was  shown,  by  mass  spectrometric 
work  that  long  lived  complex  ions  are  formed  between 
ions  and  neutral  molecules  (38) .  Recently  in  a  study  of 
the  radiolysis  of  ethanol  vapour,  Bansal  and  Freeman 
(39)  have  attributed  the  decrease  of  acetaldehyde  and 
hydrogen  yields  with  increasing  dose  to  the  scavenging 
of  the  electrons  by  acetaldehyde  with  ultimate  production 

of  CH-CHO"  (C~HrOH)  .  The  same  authors  have  postulated 

3  2  5  n 

the  formation  of  E+(CoH[.0H)  and  NH.+  (C0H_OH)  . 

A  o  n  4  a  d  n 

Formation  of  ion  clusters  in  the  D-region  (50  -  85  km 


. 
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altitude)  of  the  ionosphere  has  been  confirmed  by  the 

investigations  of  Narcisi  and  Bailey  (40) .  In  experiemnts 
carried  out  by  rocketborne  quadrupole  mass  spectrometers, 
they  have  observed  that  the  clusters  H^O+ (H^O)^  n=0 , 1 , 2 
appeared  to  be  the  dominant  positive  ions  at  an  altitude 
of  60  to  80  km.  The  "D  region  simulation"  flowing  after¬ 
glow  experiment  by  Ferguson  et  al  (41)  have  led  them  to 
suggest  a  mechanism  for  the  formation  of  the  clusters 
H^O^l^O)^  in  the  ionosphere  with  n  as  large  as  6. 

Recently,  Kebarle  et  al .  (27)  measured  the  rate  constants 

for  the  system  up  to  n=4 .  Utilizing  these 

rate  constants,  Ferguson  et  al.  (41)  calculated  con¬ 
centrations  of  Ho0+(Ho0)  near  the  concentration  maxi- 
mum  at  80  km.  They  found  that  calculated  values  agree 
well  with  the  observed  values  of  Narcisi  and  Bailey  (40) 
mentioned  above.  Thus  utilizing  Kebarle  et  al ' s  data 
the  qualitative  water  cluster  ion  production  scheme  pro¬ 
posed  earlier  has  been  made  nearly  quantitative. 

Primary  negative  ions  also  may  form  clusters  with 
neutral  molecules  in  the  D  region.  Formation  of  clusters 
such  as  °2_  ^H2°) n '  N°-(H2°)  and  C04”(H20)  have  been 
postulated  (42).  Bates  and  Dalgarno  (43)  and  W. 

Swider  (44)  have  calculated  loss  rates  for  metallic  ions 
for  various  processes  in  the  ionosphere  (9Q-100  Km  altitude) . 
They  have  postulated  that  the  lifetime  of  a  metallic  ion 
M+ (M+  =  Na,  Ca,  Mo,  K)  in  the  upper  atmosphere  above  100 
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km  is  limited  by  the  radiative  association  processes, 

M+  +  0  ->  M0+  +  hv 
M+  +  02  M02+  +  hv 

W.  Swider  (44)  suggested  that  three  body  reactions  probably 
dominate  the  loss  of  metallic  ions  near  90  km  and  below. 

+  0  +  0  2  MO  2  t  0  2 

M+  +  0  +  O2  M0+  +  C>2 

The  importance  of  alkali  ions  in  the  gas  phase 
chemistry  and  the  applications  to  the  upper  atmosphere 
are  evident  from  the  number  of  papers  devoted  to  it  at 
a  recent  symposium  (45).  The  chemistry  of  these  alkali 
ions  in  the  upper  atmosphere  is  not  yet  established 
however.  In  many  gas  phase  situations,  including  the 
upper  atmosphere  the  alkali  ions  are  the  final  ions, 
since  they  have  the  lowest  ionization  potentials  known. 
Once  the  charge  is  transferred  to  them  it  is  retained. 

We  would  suggest  that  if  the  gas  phase  contains  traces  of 
water  (as  in  the  D-region  where  H+ (Hp)  and  H  (up)  2 
account  for  35%  of  the  total  positive  ions),  the  alkali 
ions  M+  will  probably  form  MpHpp.  Like  alkali  ions, 
which  are  the  final  ions  in  the  ionosphere  due  to  their 
lowest  ionization  potential,  Ferguson  et  al  (46) 
postulated  that  N03“  ions  are  final  negative  ions  in  the 
ionosphere  because  they  are  believed  to  have  the  highest 
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electron  affinity  (46).  Recent  rocket  measurements  of 
negative  ions  in  the  ionosphere  suggest  that  NC>3  (H20)n 
(n=0-5)  may  be  the  dominating  negative  charge  carriers  (47). 
In  view  of  the  importance  of  NO^  ^H2°^n  the  i°nosP^ere ' 
a  study  is  currently  being  conducted  in  this  laboratory 
(48)  to  determine  the  thermodynamics  of  the  reactions 

N°3_ (H20) n_i  +  H2°  ?  N03~(H20)n 
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2.  EXPERIMENTAL 

2.1  General 

Measurements  were  made  with  a  60° ,  15  cm  radius  mass 
spectrometer,  which  had  been  specially  modified  for  the 
study  of  ion  hydration  in  the  gas  phase.  The  primary  alkali 
ions  were  generated  by  thermal  ionization  from  a  platinum 
filament,  coated  with  a  suitable  salt  containing  the  desired 
ion.  The  ions  react  with  water  vapour  in  the  electrically 
field  free  high  pressure  reaction  chamber  (pressures  up  to 
3  torr) .  The  residence  time  of  the  ions  is  sufficiently 
long  for  the  achievement  of  equilibrium  with  the  surrounding 
water  vapour.  Some  of  the  ions  escape  from  the  reaction 
chamber  by  diffusion  and  mass  flow  into  a  high  vacuum 
region  containing  the  accelerating  electrodes.  High  vacuum 
(~  10”4  torr)  was  maintained  by  the  low  conductance  of  the 
gas  out  of  the  reaction  chamber  and  by  the  high  capacity 
pumping  system  (pumping  speed  approximately  350  1/sec 
outside  the  chamber) . 

The  ions  were  accelerated  and  focussed  by  a  modified 
Nier  type  acceleration  system.  Detection,  after  magnetic 
mass  analysis  was  obtained  by  a  17  stage  electron  multiplier 
and  a  Victoreen  femtometer.  The  mass  spectrometer  was 
designed  and  assembled  by  G.  Collins  (49)  and  the  apparatus 
was  modified  by  S .  K.  Searles  for  the  study  of  the  gas  phase 
hydration  of  potassium  ions  (50) .  Features  of  the  instru¬ 
ment  of  particular  importance  to  the  gas  phase  hydration 
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work  are  described  in  the  subsequent  experimental  section. 

2 . 2  Instrumental  Requirements  to  Observe  Gas  Phase 

Hydration  Equilibria  and  the  Principles  of  Method. 

The  ins trumentation  required  for  the  study  of  the 
gas  phase  hydration  of  ions,  must  produce  ions,  allow  them 
to  react  and  reach  equilibrium  at  high  water  pressures  (in 
the  torr  range)  and  mass  analyze  the  resulting  ions.  To 
satisfy  these  requirements  certain  modifications  had  to  be 
made  to  a  conventional  mass  spectrometer.  A  diagram  of  the 
main  mass  spectrometer  set— up  used  in  the  present  study  is 
shown  in  Figure  2.1.  A  more  detailed  diagram  of  the  ion 
source  is  shown  in  Figure  2.2  and  2.3.  Theory  of  the  mass 
spectrometer  is  not  new  and  may  be  found  readily  elsewhere 

(51)  . 

Most  of  the  modifications  apply  to  the  reaction 
chamber  (ion  source)  of  the  mass  spectrometer.  In  order 
to  assure  the  achievement  of  equilibrium,  the  reaction 
chamber  has  to  be  able  to  handle  pressures  in  the  torr 
range  while  maintaining  the  high  vacuum  conditions  in  the 
analyzing  section.  Since  the  primary  ions  are  generated 
by  thermionic  emission  and  the  thermionic  emission  can  be 
obtained  only  at  elevated  temperatures,  the  thermionic 
source  must  be  separated  from  the  reaction  chamber  whose 
temperature  is  to  be  controlled  independently.  See  Figure 
2.2.  It  is  also  necessary  that  the  attachment  reactions 
be  carried  out  in  a  field  free  region.  Therefore,  the 
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Secondary  Electron  Multiplier  (SEM) 


FIGURE  2.1  CROSS-SECTIONAL  VIEW  OF  MASS  SPECTROMETER 
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repeller  voltage,  which  is  an  essential  part  of  the  con¬ 
ventional  mass  spectrometer  has  to  be  eliminated  from  the 
reaction  chamber  of  the  mass  spectrometer  used  in  the 
solvation  studies.  In  order  that  the  high  pressure  of 
the  reaction  chamber  does  not  interfere  with  the  performance 
of  the  mass  spectrometer,  a  high  capacity  pumping  speed 
system  is  required.  The  gas  was  pumped  through  a  six  inch 
liquid  nitrogen  baffle  and  six  inch  gate  valve  by  a 
1400  1/sec  diffusion  pump,  backed  by  a  425  1/minute  fore 
pump.  These  components  were  all  mounted  directly  below  a 
six  inch  diameter  tube  which  surrounded  the  mass  spectro¬ 
meter  reaction  chamber.  See  Figure  2.1.  The  pumping  speed 
at  the  reaction  chamber  was  measured  as  350  1/sec  for 
ammonia  with  the  liquid  nitrogen  baffle  uncooled  (50) .  The 
speed  was  calculated  from  the  pressure  drop  in  the  gas 
handling  plant  reservoir  with  respect  to  time  and  from 
the  pressure  in  the  vacuum  housing  (measured  by  an  ion 
gauge) .  In  this  measurement  the  ion  gauge  pressure  was 
2.1  x  10-4  at  a  reaction  chamber  pressure  of  3  torr .  Cooling 
the  baffle  with  liquid  nitrogen  decreased  the  ion  gauge 
pressure  to  3.8  x  10  ^  torr .  This  indicates  a  pumping 
speed  increase  to  1900  1/sec  due  to  condensation  of  ammonia 
on  the  baffle.  A  more  detailed  description  of  the  high 
capacity  pumping  system  was  given  in  references  (49,50) . 
Figure  2.2  and  2.3  show  the  design  of  thermionic  source 
used  in  the  present  work  and  the  previous  investigation 


Reaction  chamber  milled  out  of  solid  stainless 
block 

Filament 

Filament  housing 
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Cu  heat  shield  (silver  soldered  to  ion  source) 
Ceramic  insulator 
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a)  vapor  supply 
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Ion  exit  slit  leading  to  ion  accelerating  electrodes 
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FIGURE  2.2  Reaction  chamber  and  acceleration  plates  for 

hydration  reactions  of  alkali  ions 
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for  hydration  of  the  potassium  ion  (50) .  The  design  of  two 
additional  ion  sources,  which  were  necessary  to  prove 
that  the  reactions  were  not  in  thermal  disequilibrium  due 
to  the  heating  of  the  reaction  chamber  by  the  filament 
and  the  presence  of  the  drift  field  are  described  by 
S.  K  Searles  in  his  Ph.D.  thesis.  The  Kunsman  source 
(52) ,  described  in  the  next  section,  was  found  to  be  the 
most  satisfactory  source  for  all  alkali  ions  except 
lithium.  The  thermal  ionization  in  the  Kunsman  source 
arises  from  a  specially  prepared  filament  located  in  a 
housing  which  was  attached  to  the  reaction  chamber.  See 
Figure  2.2. 

The  heater  current  to  the  filament  in  Figure  2.2 
was  controlled  by  two  variable  autotransformers  in  series 
and  regulated  by  a  constant  voltage  transformer.  The 
constant  voltage  transformer  received  110  V.A.C.  from 
a  3000  V  isolation  transformer,  which  was  plugged  into 
mains,  110  V.A.C.  The  power  required  to  give  a  sufficient 
emission  of  ions  from  the  filament  at  water  vapour  pressures 
up  to  1  torr  was  about  20  W.  The  distance  between  the 
filament  and  the  reaction  chamber  was  2  cm.  A  ceramic 
insulator  thermally  and  electrically  isolated  the  filament 
housing  from  the  reaction  chamber.  Thermal  radiation  from 
the  filament  was  reflected  by  a  nichrome  reflector  backed 
by  a  copper  plate  which  dissipated  heat  energy  to  the 
reaction  chamber.  Because  the  ions  which  were  emitted  by 
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the  heated  filament  diffuse  rapidly  to  the  wall,  it  was 
necessary  to  attract  the  ions  to  the  reaction  chamber  by 
potential  difference  to  obtain  a  detectable  ion  current 
after  mass  analysis.  Experimentally,  the  highest  M 
currents  to  the  reaction  chamber  occurred  with  the  filament 
housing,  and  the  filament  at  the  same  potential.  A 
potential  difference  of  225  V  between  the  filament  housing 
and  the  reaction  chamber  was  necessary  to  attract  a  suf¬ 
ficient  number  of  ions  into  the  reaction  chamber. 

This  potential  difference  gives  an  electric  field 
strength  of  about  110  V/cm. 

2 . 3  The  Reaction  Chamber 

Figure  2.2  shows  the  reaction  chamber  attached  to  the 
filament  housing.  The  purpose  of  the  reaction  chamber  was 
to  provide  a  high  pressure  system  where  clustering  equil¬ 
ibrium  conditions  could  be  achieved  at  uniform  temperatures 
which  could  be  varied  over  a  wide  range. 

The  required  reaction  chamber  pressure  for  investi¬ 
gation  of  gas  phase  solvation  reactions  can  be  approxim¬ 
ately  calculated.  To  obtain  clusters  as  large  as  six 
ligands  about  the  central  ion,  the  ion  must  undergo  six 
collisions  plus  additional  collisions  to  remove  the  excess 
energy  produced  by  the  attachment  reaction.  For  a  dis¬ 
tance  1  of  0.5  cm  (the  distance  between  the  ion  entrance  and 
exit  slit  of  the  reaction  chamber) ,  an  ion  in  a  static  gas 
undergoes  approximately  150  collisions  at  0.1  torr  (50). 
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Since  the  ion  undergoes  approximately  15Q  collisions  at 
0.1  torr,  the  ionic  clusters  are  probably  in  equilibrium. 

To  test  the  equilibrium  assumption  it  is  desirable  to  be 
able  to  vary  the  pressure  in  the  reaction  chamber  by 
about  a  factor  of  ten,  that  is  to  increase  the  reaction 
chamber  pressure  to  1  torr. 

However,  the  mass  spectrometer,  ion  acceleration  and 
mass  analysis  region  must  be  maintained  at  a  pressure  less 
than  10  _4  torr.  Above  this  pressure  the  ion  beam  coming 
out  from  the  reaction  chamber  is  seriously  scattered  by 
collisions.  Therefore,  a  pressure  differential  of  about 
104  must  be  maintained  between  the  reaction  chamber  and 
the  ion  acceleration  region  outside.  This  can  be  achieved 
by  high  capacity  pumping  discussed  in  the  previous 
section  and  by  reducing  the  leakage  of  gas  out  of  the 
reaction  chamber.  In  the  present  instrument  the  total 
area  of  the  reaction  chamber  slots  have  been  reduced  by 
about  a  factor  of  100  and  the  pumping  speed  at  the  reaction 
chamber  increased  by  about  100  over  a  conventional  mass 
spectrometer.  Because  the  conductance  of  a  slit  increases 
when  the  mean  free  path  of  the  gas  is  smaller  than  the 
slit  width  it  was  desired  to  make  the  slit  width  less  than 
30  y,  the  mean  free  path  of  water  at  1  torr  and  273°K  (53) . 
The  increase  in  conductance  is  due  to  a  transition  from 
molecular  flow  to  viscous  flow  (54) .  The  area  of  the  ion 
exit  slit  was  Q.012  mm  x  3  mm,  typically.  The  ion  exit 
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Cross-section,  top  and  bottom  view  of  reaction 
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slit  was  made  by  spot  welding  stainless  steel  razor  blade 
edges  to  a  plate  which  was  1.5  mm  thick  and  contained  a 
1  x  4mm  centre  slot.  A  plate  on  which  the  ion  exit  slit 
was  mounted  was  bolted  to  the  bottom  of  the  reaction 
chamber  by  two  stainless  steel  screws.  See  Figure  2.3. 

The  reaction  chamber  (see  Figure  2.3)  was  machined 
from  a  single  stainless  steel  block  in  order  to  obtain 
a  more  uniform  temperature.  The  heater  holes  were  placed 
1.5  cm  from  the  reaction  chamber  cavity.  The  reaction 
chamber  lid  was  carefully  machined  to  give  a  good  fit  to 
the  reaction  chamber  bottom.  The  contact  between  the 
reaction  chamber  top  and  bottom  was  improved  by  securely 
fastening  the  two  parts  together  with  four  heavy  bolts. 

2 . 4  The  Reaction  Chamber  Pressure  Measurement 

The  pressure  of  the  reaction  chamber  was  measured  by 
means  of  a  dibutyl  phthalate  U  type  manometer  connected  with 
a  separate  line  to  the  reaction  chamber.  At  room  temperature 
the  vapor  pressure  of  dibutyl  phthalate  is  almost  negligible f 
(1mm  at  148°C) .  To  obtain  accurate  pressure  measurements 
from  the  manometer,  a  cathetometer  was  used  to  read  the  height 
difference  to  +  0.1  mm  which  corresponds  to  +  0.0077  torr 
error  in  the  ion  source  pressure. 

2 . 5  The  Reaction  Chamber  Temperature  Control  and  Measure¬ 
ment  . 

The  reaction  chamber  temperatures  above  room  temperature 
were  achieved  by  two  heaters  located  in  the  heater  holes 
drilled  in  the  wall  of  the  reaction  chamber  (See  Figure  2.3). 
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Each  heater  was  made  by  winding  an  Q.010  inch  diameter  Mo 
wire  around  a  quartz  rod  and  placing  the  rod  and  wire  coil 
inside  a  5  mm  diameter  quartz  tube.  About  3QW  was  needed 
to  heat  the  ion  source  to  40Q°C.  For  reaction  chamber 
temperatures  less  than  room  temperature,  the  heaters  were 
removed  and  the  heater  holes  were  used  as  the  cooling 
channels  of  a  circulatory  system.  Two  Kovar  seals  were 
soldered  to  the  two  heater  holes  on  the  same  face  of  the 
reaction  chamber.  Each  of  these  seals  were  attached  to 
glass  tubing  which  led  out  of  the  vacuum  housing  through 
"0"  ring  seals.  To  complete  the  reaction  chamber  cir¬ 
culatory  system  copper  tubing  was  connected  to  the  remaining 
open  ends  of  the  heater  holes.  Either  water  or  acetone,  de¬ 
pending  on  the  desired  final  temperature,  was  circulated 
through  the  reaction  chamber  by  a  Colora  Ultra-Thermostat 
which  regulated  the  bath  temperature  to  +  0,01°C.  To 
achieve  temperatures  below  0°C,  the  bath  was  filled  with 
acetone  and  cooled  with  a  2L  stainless  steel  beaker  which 
was  filled  with  dry  ice.  The  greatest  depth  of  beaker 
immersion  in  the  bath  allowed  temperatures  of  the  reaction 
chamber  as  low  as  -50 °C. 

The  temperature  of  the  reaction  chamber  was  measured 
by  means  of  an  iron-constantan  thermocouple.  The  thermo¬ 
couple  was  placed  in  a  well  drilled  through  to  the 
reaction  chamber  cavity.  The  cavity  was  capped  by  a  piece 
of  0.005  inch  thick  nichrome  sheet  which  was  spot  welded 
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to  the  wall. 

The  thermal  emf  was  measured  relative  to  a  0°C 
reference  thermocouple  by  a  Honeywell  potentiometer. 

The  temperature  of  the  thermocouple  could  be  read  to  within 
1°C. 

2 . 6  Creation  of  Alkali  Ions  and  Manufacture  of  a 
Thermionic  Source. 


It  is  well  known  that  an  abundant  supply  of  alkali 
ions  can  be  obtained  by  thermionic  emission  from  a  heated 
filament  coated  with  suitable  salt  containing  the  desired 
ion.  The  relative  probability  of  emission  of  positive  ions 
is  given  by  the  Langmuir -Saha  equation  (52) 


n 


n 


+  ■  Q -I 

a  exp  — - 
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where  n+  is  the  number  of  positive  ions  M  and  nQ  is  the 
number  of  neutral  species  M.  (j>  is  the  work  function  of  the 
filament  surface  and  I  is  the  ionization  potential  of  M. 
Because  I— 1|>  >>  kT  at  300  °K,  thermal  ionization  sources 
must  be  operated  at  high  temperatures.  The  difficulty  in 
using  a  high  temperature  source  of  ions  is  that  the  source 
has  to  be  sufficiently  distant  from  the  reaction  chamber 
so  that  the  reaction  chamber  is  relatively  unaffected.  (It 
was  desired  to  study  hydration  equilibrium  at  temperatures 
as  low  as  260°K).  The  Kunsman  source,  described  in  the 
following  paragraph  was  found  to  be  the  most  satisfactory 
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source  for  all  alkali  ions  except  lithium. 

The  manufacture  of  a  Kunsman  thermionic  source  is 
described  in  a  review  by  Jenkins  and  Trodden  (52) .  In  this 

laboratory  a  Kunsman  source  was  first  utilized  by  S .  K. 
Searles  who  investigated  hydration  of  potassium  ion  (50) . 

The  procedure  followed  in  the  present  experiment  was  as 
follows.  A  mixture  of  9  8%  Fe^,  1%  Al^  and  1%  MN03 , 
where  MN03  represents  the  corresponding  alkali  nitrate 
(M  =  Na,  K,  Rb ,  Cs)  was  fused  in  a  crucible.  In  the  fused 
product  the  Fe203  becomes  oxidized  to  Fe^.  The  product  was 
then  ground  to  a  powder.  A  Pt  filament  was  made  from  a 
strip  of  Pt  1  mm  wide  and  2  cm  long,  which  was  cut  from  a 
0.001  inch  thick  Pt  foil.  The  Pt  strip  was  spot  welded  to 
two  nichrome  studs  which  were  spaced  about  1  cm  apart.  Then 
the  ground  powder  was  fused  onto  the  filament.  This 
required  special  care  because  the  melting  point  of  the 
powder  was  close  to  the  melting  point  of  Pt.  A  thin 
layer  of  powder,  approximately  1  mm  high  was  first  spread 
uniformly  along ^he  length  of  the  filament.  A  torch  pro¬ 
vided  with  a  3/4  mm  hole  tip  was  adjusted  to  give  a  fine 
intense  flame.  The  powder  was  heated  until  it  just  began 
to  melt  and  flow.  Cycles  of  melting  and  cooling  were  used 
to  obtain  a  uniform  fused  layer  and  prevent  the  Pt  from 
melting.  If  the  coating  is  not  uniform  over  the  filament 
surface,  the  surface  temperature  and  ion  emission  will  be 
uneven.  The  reduction  of  the  surface  was  the  next  step.  The 
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chamber  was  evacuated  and  filled  with  3-10  torr  .  An 
electric  current  was  passed  through  the  filament  until  a 
dull  red  light  was  just  visible.  The  necessary  wattage  was 
15-30  W  depending  on  the  pressure.  The  emitted  ion 
current  was  collected  on  a  plate  with  a  negative  bias. 
Increasing  the  wattage  until  the  filament  became  dull  red 
yielded  a  positive  current  of  1-2  yA.  If  the  filament 
temperature  was  increased  to  a  bright  red,  a  maximum  of 
10  yA  could  be  obtained.  This  high  current  however  lasted 
only  about  20  minutes  after  which  the  current  dropped  well 
below  1  yA.  The  surface  of  the  filament  was  sufficiently 
reduced  in  about  30  minutes.  A  longer  reduction  period  did 
not  improve  the  filament  emission.  A  satisfactory  Li+  source 
could  not  be  prepared  by  the  above  procedure.  Therefore 
a  literature  search  for  suitable  thermal  ionization  source 
of  Li+  ions  was  made.  It  was  decided  to  use  3-eucryptite 
(lLi^O.  lAl^O^ .  2SiC>2 )  as  suggested  by  Blewett  and  Jones 
(55).  3-eucryptite  was  made  by  melting  Li^CO^,  AlCNO^)^  9H2° 
and  pure  powdered  quartz  together  at  about  1400 °C.  The 
melts  were  then  powdered  and  remelted  twice.  The  compound 
was  finally  powdered  and  mounted  as  a  thin  layer  on  a 
1  cm  x  0.3  cm  filament  made  of  platinum  gauze. 

2.7  Production  and  Transport  of  the  Hydrated  Ions  Out  of 

the  Reaction  Chamber. 

The  alkali  ions  are  created  in  the  filament  chamber 
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which  is  at  the  same  pressure  as  the  reaction  chamber. 

After  having  been  ejected  from  the  heated  filament,  the 
alkali  ions  drifted  to  the  reaction  chamber  under  the 
influence  of  electric  field  of  approximately  100  V/cm. 

Lower  fields  gave  insufficient  intensity.  Upon  entering 
the  field  free  reaction  chamber  the  alkali  ions  have  excess 
energy  due  to  the  drift  field.  In  the  previous  work  on 
hydration  of  potassium  ions  it  has  been  calculated  that  K+ 
ions  entering  the  reaction  chamber  possess  approximately 
0.74  eV  translational  energy.  Since  the  filament  housing 
is  at  the  same  pressure  as  the  reaction  chamber  the  alkali 
ions  are  already  partly  hydrated  when  they  enter  the  reaction 
chamber.  Thus  in  spite  of  relatively  high  potential  dif¬ 
ference  between  the  filament  and  the  reaction  chamber 
the  ions  entering  the  reaction  chamber  possess  relatively 
low  (~  0.74  eV)  translational  energy.  In  the  reaction 
chamber  the  ions  are  further  thermalized  by  losing  the 
excess  translational  energy  through  collisions  with 
solvating  molecules.  After  the  primary  ions  lose  their 
excess  energy,  the  ions  react  with  neutral  molecules  as  they 
diffuse  toward  the  wall.  During  this  time  the  ions  should 
form  equilibrium  ionic  clusters.  On  their  way  towards  the 
walls  of  the  reaction  chamber,  where  they  are  discharged, 
the  primary  ions  collide  with  the  solvating  molecules.  Some 
of  these  collisions  may  result  in  the  attachment  of  the 

to  the  neutral  molecules  thus  forming  ion-clus ters .  The 
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ion-clusters  stabilize  themselves  by  giving  up  their  excess 
energy  through  collisions  with  other  molecules  (third  body 
collisions) .  In  order  that  equilibrium  conditions  prevail, 
this  stabilization  of  ion  clusters  should  be  achieved, 
before  the  ion  clusters  diffuse  out  of  the  reaction  chamber. 

The  time,  t,  necessary  for  the  ions  to  travel  the 
distance  d  can  be  estimated  by  the  approximate  equation  (I) 
given  below,  which  establishes  a  relationship  between  the 
time,  t,  required  by  a  particle  to  diffuse  the  distance,  d 
(56),  and  the  diffusion  coefficient,  D. 


2D 


In  the  present  case  d  is  the  distance  between  the  ion 

entrance  and  exit  slit  of  the  reaction  chamber  and  is 

equal  to  1.2  cm,  D  is  the  diffusion  coefficient  of  the 

2 

ions  through  the  gas  and  is  taken  as  500  cm  /sec  (56)  for  an 
ion  source  pressure  of  0.1  torr  and  temperature  300 °K. 

Thus,  the  diffusion  time  t  is  calculated  to  be  2.9  milli- 
sec.  This  time  is  an  estimate  of  the  ion  residence  time 
in  the  reaction  chamber.  This  estimate  is  approximate 
since  space  charge  effects  due  to  presence  of  the  charged 
particles  may  alter  the  diffusion  coefficient. 

The  half-life  of  the  attachment  reaction 

2.1  M+  +  2H20  +  M+ (H20)  +  H20 
can  be  calculated  if  we  assume  a  value  for  the  rate  constant 


k.  The  reaction 
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2.2  Cl  +  2H20  +  Cl  (H20)  +  H20 

has  been  studied  by  Payzant  and  Kebarle  (57)  and  Lineberger 

and  Puckett  (58) .  They  found  the  third  order  rate  con- 

-29  2  2  . 

stants  to  be  k  ;  10  cc  /molecule  sec.  As  an  approximation 

the  same  value  can  be  assumed  for  the  rate  constant  for  the 

reaction  2.1,  because  the  ion  M+  is  monoatomic  and  reaction 

2.1  bears  resemblance  to  reaction  2.2.  The  half-life  t, 

£ 

is  then  given  by 


II 


1 

TTq  2 

10  [H20] 


where  [H20]  is  the  concentration  of  water  in  the  reaction 

chamber  in  molecules/cc.  At  a  reaction  chamber  pressure 

of  0.1  torr  and  temperature  300 °K,  the  water  concentration 

is  equal  to  0.3  x  10  molecules/cc.  Introducing  this 

value  into  equation  (II) ,  the  half-life  is  found  to  be 

~1  milli-sec.  This  is  roughly  a  factor  of  three  smaller 

than  the  ion  residence  time  calculated  from  equation  (I) . 

d2 

Thus  the  reaction  time  (the  ion  residence  time  t  ~  ) 

would  appear  to  be  sufficiently  long  even  at  0.1  torr  for 
the  reaction  2.1  to  reach  an  equilibrium.  However  in 
the  present  studies  the  lowest  reaction  chamber  pressures 
used  were  0.3  torr. 

In  addition  to  the  diffusion  process,  the  ions  can 
be  removed  from  the  ion  source  by  mass  flow.  It  has  been 
shown  (50) ,  however,  that  the  mass  flow  does  not  become 
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significant  until  about  1  mm  from  the  exit  slit.  Thus 
starting  from  the  entrance  to  the  reaction  chamber,  the 
ions  are  transported  for  the  first  3  to  4  mm  by  diffusion 
and  in  the  last  1  mm  by  mass  flow. 

2 . 8  Ion  Acceleration,  Mass  Analysis  and  Ion  Current 

Detection . 

Ions  diffusing  out  of  the  ion  exit  slit  of  the 
reaction  chamber  were  accelerated  from  the  reaction  chamber 
potential,  +  1850  V,  to  earth  potential.  Figure  2.2  shows 
the  electrode  system,  which  was  a  modified  Nier  design. 

The  first  electrode  spacing  to  the  reaction  chamber  was 
increased  from  1  mm  to  10  mm.  As  discussed  in  detail  in 
reference  (50),  the  electrode  when  present  with  the 
shorter  distance  presented  a  pumping  obstacle  to  the  gas 
escaping  from  the  ion  exit  slit.  The  resulting  high 
pressure  caused  the  ionic  clusters  to  undergo  collisions 
outside  the  reaction  chamber,  which  led  to  non- eguilibr ium 
effects.  After  acceleration  the  ions  entered  the  magnetic 
analysis  tube  and  were  deflected  by  60°.  The  resolution 
calculated  from  the  magnetic  analyzer  entrance  and  exit 
slit  widths  (49)  was  200.  The  observed  resolution  was 
slightly  less. 

After  the  ions  passed  through  the  analyzer  exit  slit, 
they  entered  the  17  stage  secondary  electron  multiplier 
(SEM) .  In  the  SEM  a  positive  ion  is  further  accelerated 
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by  the  potential  of  the  first  dynode.  On  striking  the 
dynode ,  secondary  electrons  are  ejected  and  accelerated 
to  the  second  dynode  where  each  electron  ejects  more 
electrons  etc.  The  current  amplification  of  this  device 
is  a  function  of  the  voltage  across  its  dynodes .  Normally 
there  is  a  uniform  potential  drop  between  each  dynode,  the 
last  dynode  (collector)  being  at  ground  potential.  The 
ratio  of  the  electron  current  from  the  17th  dynode  of  the 
SEM  to  the  positive  ion  current  collected  from  the  first 
dynode  is  the  absolute  gain  of  the  SEM.  The  currents 
were  measured  by  a  Victoreen  Femtometer  with  a  10^  ohm 
input  resistor.  In  the  case  of  the  M+tH^O)  series  the 
ion  currents  were  too  small  to  be  measured  on  the  first 
dynode,  therefore  the  absolute  gain  for  M^K^O)  series 
could  not  be  determined. 

2.9  Relation  of  the  Ion  Current  to  the  Ion  Concentration 
From  the  time  of  their  formation  in  the  reaction 
chamber  to  the  time  of  their  recording  by  the  detection 
system  the  ion-clusters  are  engaged  in  a  number  of  pro¬ 
cesses.  A  few  of  these  processes  may  discriminate  against 
one  or  the  other  type  of  ions.  Some  of  these  discrimin¬ 
atory  processes  will  be  discussed  below. 

A.  Discrimination  in  transport  through  the  reaction  chamber 
The  diffusion  coefficient  of  an  ion  is  inversely  pro¬ 
portional  to  the  square  root  of  the  mass  of  the  ion  (56) . 
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Therefore  one  expects  that  the  lighter  ions  would  travel 
through  the  reaction  chamber  faster  and  would  reach  the 
ion  exit  slit  in  a  shorter  time  than  the  heavier  ions. 

Thus,  the  actual  concentration  of  the  ions  with  smaller 
mass  would  be  less  than  what  is  represented  by  the  recorded 
signal.  However,  since  the  average  half-life  of  the 
clustering  reaction  is  short  in  comparison  to  the  residence 
time  of  the  ion  and  thus  a  cluster  M+ (H2°) n  gains  and 
loses  water  molecules  at  a  fast  rate;  all  the  clusters  would 
have  the  same  average  molecular  weight  and  their  diffusion 
coefficient  would  be  equal.  However  at  the  moment  when 
the  ions  effuss  through  the  leak  Cion  exit  slit)  ,  discrimin¬ 
ation  might  be  expected  to  occur  since  effusion  rate  through 


the  leak  is  proportional  to  the  average  thermal  velocity , 
v  oc  The  A  difference  is  highest  and  equal  to  1.88 

v  ^  24 


M 


for  Li+  and  Li+H20  ratio.  However  this  reaction  was  not 


T 

M 


for  the  Na+  and  Na  H20  is 


+ 


measured.  The  difference  in 

equal  to  1.33.  Thus  KQ  for  the  reaction:  Na"  +  H20  + 
Na+(H20)  might  be  out  by  a  factor  of  1.33.  This  correction 
is  small  and  has  been  neglected. 


B .  Discrimination  due  to  stripping  outside  of  the_ 

reaction  chamber. 

If  the  pressure  in  the  region  just  outside  of  the  ion 
exit  slit  is  high,  some  ion  molecule  collisions  may  take 
place  in  the  vacuum  region  between  the  reaction  chamber 
exit  slit  and  the  first  accelerating  electrode.  This 
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problem  was  encountered  by  Collins  (49)  in  a  preliminary 
study  of  the  equilibrium  reaction 

H+(H20)n-l  +  H2°  ?  H+(ll20)n 

Collins  found  that  the  ion  cluster  size  began  to  decrease 
after  a  pressure  of  about  0.3  torr  due  to  ion  cluster- 
water  molecule  collisions  outside  the  reaction  chamber. 

Due  to  the  presence  of  the  ion  accelerating  field,  these 
collisions  occur  with  kinetic  energies  much  higher  than 
thermal.  As  a  result  of  these  collisions,  the  ion- 
clusters  may  dissociate.  Since  a  larger  cluster  has 
weaker  ion-ligand  bond  strengths  and  may  have  larger 
collision  cross  section,  the  larger  clusters  were  prefer¬ 
entially  dissociated.  The  extraneous  collisions  were 
eliminated  by  increasing  the  pumping  speed  at  the  ion  exit 
slit  by  increasing  the  distance  between  the  reaction  chamber 

and  the  first  electrode. 

As  shown  in  Chapter  3,  the  calculated  equilibrium 
constants  did  not  change  with  pressure  up  to  some  maximum 
pressure  P  after  which  the  equilibrium  constant  began  to 
decrease.  Evidently  at  higher  pressures  nonequilibrium 
collisions  again  became  significant. 

C #  Discrimination  in  Mass  Analyzer  Tube — 

In  general,  in  a  magnetic  mass  analyzer  tube  such  as 
the  one  used  in  this  work  any  loss  of  ions  due  to  elastic 
scattering,  collisions  with  the  wall  etc.,  will  be  the 
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same  for  all  ions.  Consequently,  we  can  assume  that  the 
effect  of  mass  discrimination  in  the  analyzer  tube  was 
negligible . 

D.  Discrimination  in  the  ion  current  detection  system 
Discrimination  may  also  occur  in  the  secondary 
electron  multiplier.  The  number  of  electrons  ejected  from 
the  first  dynode  of  SEM,  is  dependent  on  the  velocity  of 
the  ions  as  well  as  their  chemical  composition  (59) .  Thus, 
the  gain  of  the  multiplier  would  not  be  the  same  even  for 
two  isobaric  ions  such  as  N+  and  CH2  +  (60)  .  The  reason 
for  this  difference  is  that,  the  multi-atom  ions  may  break 
up  when  they  strike  the  surface  of  the  SEM  dynode  and  act 
as  a  group  of  particles  instead  of  just  one  ion.  To 
measure  the  SEM  discrimination  the  relative  gain  of  the 
SEM  for  the  series  of  solvated  ions  should  be  measured. 
However,  in  the  case  of  the  M+(H20)n  series  the  relative 
gain  of  the  SEM  could  not  be  measured,  since  the  ion  currents 
on  the  first  dynode  were  too  small  to  be  measured  directly. 

From  the  work  of  S.  K.  Searles  (5Q)  the  relative 
gains  for  H+(H20),  H+(H20)2,  H  (H20) 3  and  H  (H20) 4  were 
1.0,  1.1,  0.8,  0.8  respectively.  For  NH^  (NH^) ,  NH^  (NH^) 2 
and  NH4+(NH3) 3  the  relative  gains  were  1.0,  0.95,  0.95 
respectively.  Therefore,  it  appears  that  the  SEM  introduces 
a  maximum  discrimination  of  2Q%  in  the  measurement  of  ionic 
clusters.  This  error  is  considered  small  since  the  AG 
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depend  on  the  log  K  ,  and  are  thus  little  affected  by  a 

n— i,n  u 

20%  error.  The  AH  ,  are  completely  unaffected  since 

n-1 , n  J 

they  do  not  change  by  the  presence  of  a  constant  factor. 
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3 .  RESULTS  OF  HYDRATION  OF  POSITIVE  ALKALI  IONS  M+ 

3 . 1  Introduction  and  Thermodynamic  Relations 

Alkali  ions  which  are  ejected  by  a  heated  filament 
drift  into  the  ion  source  under  the  influence  of  an  elect¬ 
ric  field.  In  the  field  free  ion  source,  the  alkali  ions 
lose  their  excess  translational  energy  and  react  with  water 
molecules  as  follows: 

M+  +  H20  y  M+H20 

m+(h2o)  +  h2o  +  m+(h2o)2 


(H20) n-l  +  H2°  t  M+(H20)n 


Provided  that  the  equilibrium  is  achieved  in  the  ion  source 
the  equilibrium  constant  can  be  calculated  from 


K 


M+  (Ho0) 

2  n 


n-l ,  n 


PM+  (H2°)  n-l  Pfi2° 


where  the  corresponding  equilibrium  concentrations  of 

M+(u  0)  and  M+(Ho0)  ,  are  assumed  to  be  proportional  to 

2  n  2  n-l 

the  mass  analyzed  ion  currents.  Determination  of 

O 

at  different  temperatures  is  used  to  obtain  AH  n_ifUr 


AS  ,  and  AG 

n-l,n  n-l,n 
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3.2  Presentation  of  the  Results  and  General  Discussion 
The  equilibrium  constants  obtained  for  the  gas 
phase  hydration  of  the  lithium,  sodium,  rubidium  and  cesium 
ions  are  presented  in  Figures  3.1  to  3.21.  It  can  be  seen 
that  at  any  one  temperature  the  values  of  K  ,  remain 

II  JL  /  1 1 

relatively  constant  over  the  pressure  range  covered.  This 
independency  of  Kn_^  n  of  water  pressure  indicates  that 
equilibrium  conditions  in  the  ion  source  were  achieved  and 
no  significant  amounts  of  collisions  occurred  outside  of 
the  ion  source. 

Van't  Hoff  type  plots  of  the  equilibrium  constants 
are  shown  in  Figures  3.22  to  3.25.  The  equilibrium  con¬ 
stants  presented  in  these  Figures  are  2  to  K5^  for 

lithium,  Kq  g  to  K5  6  for  sodiuin'  Ko  1  to  K4 , 5  for  rubidium 

and  _  to  K~  .  for  cesium  ion  clusters.  The  Kn  ,  values 

0,1  3,4  u  /  x 

for  lithium  ions  were  not  measured  since  the  high  tempera¬ 
tures  (T  >  9 0 0 ° K)  were  required  to  observe  the  equilibrium 

constant  K~  .  The  reaction  chamber  used  in  the  present 

U  /  J- 

study  could  be  operated  at  the  maximum  temperature  of  600 °K 

while  the  Kn  equilibrium  for  the  lithium  ion  could 

0,1 

probably  be  measured  with  some  refinements  of  the  heating 
system  of  the  reaction  chamber.  The  additional  effort 
involved  was  not  considered  necessary  for  the  present  time 

O  o 

since  a  fairly  accurate  value  for  AH  ^  ^  and  AG  q  could 

O 

be  obtained  by  extrapolation  from  the  other  AH  n_q,n  and 

O 

AG  of  this  ion. 

n-1 ,  n 

The  values  for  AH  were  determined  from  the 

n  l  /  n 
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FIGURE  3 . 1  Plots  of  Log  K ^  2  f°r  the  Gas.  Phase  Hydration 


of  Li  + 


at  Various  Temperatures,  Versus  Pressure  of  H20. 
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FIGURE  3.2  Plots  of  Log  Kn  0  for  the  Gas  Phase  Hydration 
of  Li  ,  at  Various  Temperatures,  Versus  Pressure  of  H2O. 
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FIGURE  3.3  Plots  of  Log  K3  4  for  the  Gas  Phase  Hydration 
of  Li+,  at  Various  Temperatures,  Versus  Pressure  of  H20. 
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FIGURE  3.4  Plots  of  Log  K„  c  for  the  Gas  Phase  Hydration 
of  Li  ,  at  Various  Temperatures,  Versus  Pressure  of  HO. 
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FIGURE  3.5  Plots  of  Log  Kc  c  for  the  Gas  Phase  Hydration 
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of  Li  ,  at  Various  Temperatures,  Versus  Pressure  of  H^O. 
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Plots  of  Log  Kq  ^  for  the  Gas  Phase  Hydration 
Various  Temperatures,  Versus  Pressure  of  H^O. 
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FIGURE  3.7  Plots  of  Log  K-,  0  for  the  Gas  Phase  Hydration 
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of  Na  ,  at  Various  Temperatures,  Versus  Pressure  of  H^O. 
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FIGURE  3.8  Plots  of  Log  K2  3  for  the  Gas  Phase  Hydration 
of  Na+ ,  at  Various  Temperatures,  Versus  Pressure  of  H20. 


© 


W 


T'r 


@ 


© 


119°- 


146' 


166 


200° 


230° 


LOG 


57. 


10 


CO 

**> 

CM 

0 


-1.0 


0.5  1.0 

P(i’orr) 

FIGURE  3.9  Plots  of  Log  K2  3  for  the  Gas  Phase  Hydration 
of  Na+,  at  Various  Temperatures,  Versus  Pressure  of  H20 
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FIGURE  3.10  Plots  of  Log  K3  4  for  the  Gas  Phase  Hydration 
of  Na+/  at  Various  Temperatures,  Versus  Pressure  of  H20 
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FIGURE  3.11  Plots  of  Log  ^  for  the  Gas  Phase  Hydration 
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of  Na  ,  at  Various  Temperatures,  Versus  Pressure  of  H^O 
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FIGURE  3.12  Plots  of  Log  c  for  the  Gas  Phase  Hydration 
—  ,  b 

of  Na  ,  at  Various  Temperatures,  Versus  Pressure  of  HO 
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FIGURE  3.13  Plots  of  Log  KQ  ^  for  the  Gas  Phase  Hydration 
of  Rb+ ,  at  Various  Temperatures,  Versus  Pressure  of  H20 
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FIGURE  3.14  Plots  of  Log  K-,  0  for  the  Gas  Phase  Hydration 
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of  Rb  ,  at  Various  Temperatures,  Versus  Pressure  of  H^O 
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FIGURE  3.15  Plots  of  Log  K2  3  for  the  Gas  Phase  Hydration 
of  Rb+ ,  at  various  Temperatures,  Versus  Pressure  of  H20 
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FIGURE  3.16  Plots  of  Log  ^  f°r  the  Gas  Phase  Hydration 
of  Rb+ ,  at  Various  Temperatures,  Versus  Pressure  of  H20 


■ 


LOG 


65  . 


P(forr) 

FIGURE  3.17  Plots  of  Log  K„  c  for  the  Gas  Phase  Hydration 
'  *  T  q  b 

of  Rb  ,  at  Various  Temperatures,  Versus  Pressure  of  H20 
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FIGURE  3.18  Plots  Of  Log  Kn  ,  for  the  Gas  Phase  Hydration 
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of  Cs  ,  at  Various  Temperatures,  Versus  Pressure  of  H^O 
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FIGURE  3.19  Plots  of  Log  K±  2  for  the  Gas  Phase  Hydration 
of  Cs+,  at  Various  Temperatures,  Versus  Pressure  of  H2<D 
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FIGURE  3.20  Plots  of  Log  K2  3  for  the  Gas  Phase  Hydration 
of  Cst  at  Various  Temperatures,  Versus  Pressure  of  H„0 
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FIGURE  3.21  Plots  of  Log  K3  4  for  the  Gas  Phase  Hydration 
of  Cs+,  at  Various  Temperatures,  Versus  Pressure  of  h^O 
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FIGURE  3.22  Van ' t  Hoff  Type  Plots  of  the  Equilibrium 
Constants  for  the  Gas  Phase  Hydration  of  Li+ 
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Constants  for  the  Gas  Phase  Hydration  of  Na+ 
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FIGURE  3.24  Van 1 1  Hoff  Type  Plots  of  the  Equilibrium 
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Constants  for  the  Gas  Phase  Hydration  of  Rb 
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FIGURE  3.25  Van ' t  Hoff  Type  Plots  of  the  Equilibrium 
Constants  for  the  Gas  Phase  Hydration  of  Cs+ 
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slopes  of  Van't  Hoff  plots,  by  means  of  least  square  treat¬ 
ment  using  values  of  K  ,  ,  measured  at  various  pressures 

1 1  J-  f  1 1 

and  temperatures. 


o  o 

AG  n  and  AS  ,  at  298b  were  calculated 

n-l,n  n-l,n 

O 

from  equation  (1.1)  and  (1.2)  .  The  values  of  AH  n_g^n  as 
well  as  values  of  AG°  n  „  and  AS  are  presented  in 

Tables  3.1  to  3.4  and  in  Figures  3.26  and  3.27. 

O 

In  Figure  3.26  the  AH  ,  for  the  ions  are  shown 

n  i  f  n 

for  the  different  clustering  steps  (n-l,n) .  It  is  found 

0 

for  a  given  alkali  ion,  that  -AH  n_1  n  decreases  as  the 
size  of  the  cluster  (n)  increases.  Comparing  the  differ— 

O 

ent  alkali  ions  one  finds  that  the  AH  n-1  are  largest 

for  Li+  and  smallest  for  Cs+.  Both  results  are  reasonable 

and  agree  with  what  one  would  have  expected  namely,  that 

AH  would  decrease  with  increase  of  the  size  of  the 

n-1 ,  n 

ion.  The  AG  values  in  Figure  3.27  exhibit  similar 

n  — ’  -L  f  n 

trends . 

O 

The  consistency  of  the  relative  AH  n-1  values 

for  the  different  alkali  ions  with  available  thermodynamic 

data  is  checked  in  Figure  3.28  which  gives  a  plot  of: 

AH  (Cs+)  -  AH  (M+)  versus  n.  It  is  to  be  expected  that 

o,n  o,n 

for  large  n  the  above  difference  should  approach  the  differ¬ 
ence  between  the  simple  ion  heats  of  hydration:  AH^(Cs  )  - 


AH  (M+) .  Indicated  on  the  figure  are  the  values  for  these 
h 

differences  obtained  from  the  simple  ion  hydration  enthalpies 
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TABLE  3 . 1 


Experimental 

Thermodynamic 

Values  for  the  Gas 

Phase  Reactions: 

Li+ 

(Ho0)  +  Ho0 

2  n-1  2 

t  Li+(H2o)n 

(n-1 ,n) 

Reaction 

O 

-AH  , 

n-1 ,  n 

-AG°  ,  a 
n-1 ,  n 

-AS  a 

n-1 ,  n 

£ 

i — 1 

i 

G 

(kcal/mole ) 

(kcal/mole) 

(e . u. ) 

1 — 1 

o 

34. 0b 

25. 5b 

23. 0b 

1,2 

25.8 

18.9 

23.1 

2,3 

20 . 7 

13.3 

24.9 

3,4 

16.4 

7.5 

29.9 

4,5 

13.9 

4.5 

31.4 

5,6 

12.1 

2.5 

32.0 

a 

Standard  state  1  atm  and  298°K. 

^Extrapolated  values  from  plots  in  Figures  3.26  and  3.27 


' 
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TABLE  ^  2 


Experimental  Thermodynamic  Values  for  the  Gas  Phase  Reactions: 


Na+ (H„0)  ,  +  H„0 

2  n-1  2 

i  Na+(H20)n 

c 

• — 1 

1 

G 

Reaction 

O 

-AH 

n-1 ,  n 

-AG°  a 

n-1 ,  n 

-AS  ,  a 

n  - 1 ,  n 

n-1,  n 

(kcal/mole) 

(kcal/mole) 

(e  .  u . ) 

i — 1 

o 

24.0 

17.6 

21.5 

1,2 

19.8 

13.2 

22 . 2 

2,3 

15 . 8 

9.3 

21.9 

3 , 4 

13.8 

6 . 3 

25 . 0 

4,5 

12 . 3 

3.9 

28.1 

5,6 

10.7 

2.9 

26.0 

Standard  state  1  atm  and  298°K. 
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TABLE  3 . 3 


Experimental 

Thermodynamic 

Values  for  the 

Gas  Phase  Reactions: 

Rb+ 

(Ho0)  ,  +  Ho0 

2  n-1  2 

t  Rb+(H20)n 

(n-1, n) 

Reaction 

O 

-AH 

n-1 ,  n 

-AG°  a 

n-1 ,  n 

°  3 

-AS 

n-1 ,  n 

rt-1 ,  n 

(kcal/mole) 

(kcal/mole) 

(e  .  u . ) 

0,1 

15.9 

9 . 6 

21.2 

1 , 2 

13.6 

7.0 

22 . 2 

2,3 

12.2 

5.0 

24 . 0 

3,4 

11.2 

3.8 

24.8 

4,5 

10.5 

2.8 

25.7 

a 

Standard  state  1  atm  and  298°K. 


' 
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TABLE  3 . 4 


Experimental 

Thermodynamic 

Values  for  the  Gas 

Phase  Reactions: 

Cs+(Ho0}  ,  + 

2  n-1 

H2°  J  Cs+(H20)n 

(n-1 , n) 

Reaction 

O 

-AH  , 

n-1 ,  n 

-AG°  a 

n-1 ,  n 

O  rs 

-AS 

n-1 ,  n 

n-1,  n 

(kcal/mole) 

(kcal/mole) 

(e  .  u . ) 

i — 1 

o 

13.7 

7.9 

19.4 

1,2 

12 . 5 

5.9 

22.2 

2,3 

11.2 

4.2 

23 . 7 

3,4 

10.6 

3.0 

25.4 

a 

Standard  state  1  atm  and  298°K. 
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AH  Hydration 

FIGURE  3.29  Plot  of  enthalpies  for  clustering  reactions 

M  (H2°^n-1  +  H2^  M  (I^O)  versus  total  single  ion  heats 
of  hydration  obtained  by  Randles. 


83. 


of  Latimer,  Pitzer  and  Slanski  (32)  and  Randles  (35) .  The 

plot  clearly  shows  that  the  present  data  are  consistent  with 

the  simple  ion  hydration  differences.  In  Figure  3.29 

AH°  n  values  were  plotted  versus  AH  hydration  of  Randles 
n-1 , n  c 

®  ,  -{-  -f* 

(35) .  This  Figure  shows  that  AH  for  Li  ,  Na  ,  K  , 

X 1  -L  ^  11 

Rb+  and  Cs+  ion  clusters  fall  on  the  same  curve  and  that 
0 

AH  decreases  as  the  size  of  the  ion  and  size  of  the 

n-1 ,  n 

cluster  increases. 

3.3  Equilibrium  Distribution  of  M+ (H^O) 

The  equilibrium  constant  values,  obtained  in  this 
work  can  be  utilized  to  calculate  the  relative  equilibrium 
concentrations  of  the  cluster  M  (H2O)  at  various  temper¬ 
atures  and  partial  pressures  of  water.  As  an  illustration, 
the  relative  equilibrium  concentrations  of  the  clusters 

Na+(H20)n  are  presented  in  Figure  3.30  for  a  temperature  of 

-5 

300 °K  and  partial  pressures  of  water  ranging  from  10  to  10 
torr.  For  this  calculation  the  equilibrium  constants  Kn-1  n 
for  n  =  1  to  6  were  taken  from  Figure  3 . 2  2  at  a  temperature 
of  300 °K.  The  concentration  of  Na+  was  arbitrarily  set 
to  unity.  The  concentration  of  the  cluster  Na+(H20)n  is 
given  by 

[Na+(H20)1  =  Ko<nPn 

where  K  is  the  product  of  the  equilibrium  constants 
o ,  n 

from  Kn  ,  to  K  ,  and  P  is  the  partial  pressure  of 

0,1  n-l , n 


»oOoe 
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"f* 

of  water.  By  setting  the  sum  of  the  clusters  Na  to 

Na+(Ho0)  to  100,  the  relative  concentration  of  each 
2  n 

cluster  is  obtained.  The  plots  in  Figure  3.30  show  that 
at  low  pressures  the  lower  hydrates  are  present,  the  dif¬ 
ference  between  the  stability  (free  energy)  of  the  hydrates 
is  large  and  the  concentrations  of  one  or  two  hydrated 
species  are  dominant.  At  higher  pressures  where  the 
higher  hydrates  become  stable,  differences  between  the 
free  energies  of  the  hydrates  are  small  and  a  number 
of  hydrated  species  are  present  simultaneously. 


( 
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4 •  ELECTROSTATIC  calculations  of  the  potential  energies  of 
THE  CLUSTERS  M+(H20)n  AND  COMPARISON  WITH  EXPERIMENTAL 

ENTHALPY  CHANGES 

In  this  section,  the  interaction  energy  between  an  ion 

and  its  neighbouring  water  molecules  is  discussed.  The 

energy  of  interaction  between  a  spherical  ion  and  neutral 

water  molecule  is  expressed  as  a  sum  of  energy  terms.  The 

energy  terms  represent  various  forces  which  are  known  to 

act  between  the  ion  and  the  surrounding  water  molecules. 

The  calculated  potential  energy  changes,  (AEq  which 

occur  when  the  water  molecules  bond  to  the  ion,  will  then 

be  compared  to  the  experimental  enthalpies  of  hydration 

Ah  of  the  ions  at  298°K. 
o,  n 

4 . 1  Interaction  Energies 

The  energy  of  interaction  between  an  inert-gas  like 

* 

ion  and  a  water  molecule  is  assumed  to  be  the  sum  of  the 
following  contributions. 

A.  Ion-permanent  dipole  forces:  The  attractive  energy 

between  the  ion  and  the  permanent  dipole  of  a  water 

2 

molecule  is  zeycos0/R  ,  where  y  is  the  dipole  moment, 

ze  the  charge  and  0  the  angle  between  the  dipole  and 

the  line  of  centers.  This  attractive  energy  term  has 

a  directional  character,  the  maximum  value  being 
2 

zey/R  (61) . 

B.  Ion-induced  dipole  forces:  The  ion  induces  a  dipole 

in  the  polarizable  water  molecule  and  an  attraction 

2  2  4 

results.  The  energy  being  z  e  a/2r  ,  where  a  is  the 


■ 


polarizability  (61). 


Dispersion  forces;  Even  when  species  have  no  dipole 
moments  there  is  an  attraction  between  them  due  to 
mutual  polarization.  The  interaction  energy  corres¬ 
ponding  to  these  dispersion  forces  leads  to  a  potential 
energy  (61 ) . 


E  =  3/2 


ala2 


r;  . 
1  /  3 


X1X2 


Il+I2 


The  I's  represent  the  ionization  energies  of  the  ion  and 

molecule,  the  a's  are  their  polarizabilities  and  R-  •  is 

1  r  J 

the  internuclear  separation  between  the  ion  and  the  ligand 
molecule  and  between  the  ligand  molecules  themselves. 


Dipole -dipole  repulsive  forces:  These  are  directional, 

3 

the  maximum  interaction  energy  being  2yAyB/R  when  the 
dipoles  are  in  line  with  each  other.  The  y's  are  dipole 
moments  of  molecules  and  R  is  the  distance  between  them. 


Ion-ligand  molecule  repulsive  forces:  The  ion-molecule 
repulsive  energies  are  interpreted  by  using  a  term  of 
the  form  AR~n ,  A  and  n  are  the  constants  and  R  is  the 
internuclear  separation  between  the  ion  and  the  ligand 
molecule.  A  value  of  12  for  n  is  assumed  in  the  well- 
known  Lennard-Jones  potential. 


■ 
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4 . 2  The  Present  Model  for  the  Electrostatic  Calculations 
The  potential  energy  E^  of  a  cluster  rela¬ 

tive  to  that  of  the  alkali  ion  and  n  water  molecules  at 
infinity  can  be  written  as 

4.1  Er  =  EDIP  +  EPOL  +  EDIS  +  RDIP  +  REL 

where  EDIP  represents  the  attraction  energy  between  the 
central  ion  and  the  permanent  dipole  of  the  ligand  molecule, 
EPOL  represents  the  attraction  energy  due  to  the  polariza¬ 
tion  of  the  ligand  molecules  by  the  central  ion,  EDIS 
represents  the  dispersion  energy  between  the  central  ion 
and  the  ligands  as  well  as  the  dispersion  energy  between 
the  various  ligands,  RDIP  represents  the  dipole-dipole 
repulsion  energy  between  the  ligand  molecules  and  finally, 

REL  represents  the  ion  ligand  repulsion  energy.  The 
electronic  repulsions  between  the  ligand  molecules  themselves 
were  neglected.  The  models  and  equations  used  in  the 
present  work  were  designed  to  give  a  good  approximation  to 
the  sum  of  the  electrostatic  energies,  given  by  equation 
4.1.  Each  term  in  equation  4.1  will  be  discussed  separately. 

A.  Ion -permanent  dipole  interaction  (EDIP) 

It  is  well  known  that  an  asymmetric  molecule, 

for  example,  possesses  a  permanent  electrical  dipole.  The 
permanent  dipole  of  a  water  molecule  is  due  to  a  permanent 
charge  shift  of  the  electron  cloud  in  the  water  molecule. 


■  ■  .  .  -  M  . 
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At  large  distances  between  the  ion  and  the  water  molecule, 
the  point  dipole  approximation  for  the  water  molecule  is 
completely  satisfactory.  However,  as  the  distance  to  the 
dipole  molecule  is  decreased  the  point  dipole  approximation 
becomes  unsatisfactory.  An  approximation  that  is  somewhat 
better  (at  close  distances  between  the  ion  and  dipole)  is 
to  represent  the  permanent  dipole  by  a  few  point  charges 
(62) •  The  placement  of  the  positive  and  negative  point 
charges  should  be  such  as  to  reproduce  the  point  dipole 
(at  large  distances) .  For  this  calculation  it  is  reason¬ 
able  to  place  the  negative  charge  on  the  oxygen  nucleus 
and  positive  charges  on  the  hydrogen  nuclei  of  a  water 
molecule.  For  a  given  charge  separation  the  size  or  the 
charges  is  determined  from  the  permanent  dipole  of  water 
1.85  Debye  and  the  bond  angle  of  105°.  The  O-H  bond  dis- 

O 

tance  was  taken  as  0.97  A  (63 ) • 

To  determine  the  energy  of  interaction  of  an  ion  and 
the  permanent  dipole,  one  simply  adds  the  coulontoic  terms 
due  to  the  interaction  of  the  ion  with  the  three  poini_ 
charges  of  the  water  molecule.  The  point  charge  model  used 
in  the  present  calculation  is  shown  in  Figure  4.1  and 
equation  4.2.  The  same  model  has  been  used  by  Doyne  and 
Caldwell  (62)  who  considered  the  H30+(H20)n  system. 

n  p  (334)  (Qion)  (Q.  ,) 

4.2  EDIP  =  l  l  - - - - 

.  ,  .  ,  R.  • 

3=1  i=1  i»] 
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H;  atom  1 


H,otom  2 


FIGURE  4.1 


M+ (H20) 


model  used  in  the  theoretical 


calculations . 
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.  th 

Q^on  is  the  ionic  charge,  ^  is  the  charge  on  1  atom 

,  th 

of  the  j  molecule  with  p  atoms  per  molecule  and  n  mole¬ 
cules  per  cluster.  R.  .  is  the  ion-atom  distance. 

1  /  3 

The  numerical  factor  (334)  in  the  above  equation  is 
adjusted  to  give  the  energy  in  kcal/mole ,  for  values  of 

O 

Q  in  atomic  units  and  R  in  A. 


B .  Ion-induced  dipole  interaction  (EPOL) 

The  ion,  carrying  a  charge  Q,  produces  an  electric 


field  of  magnitude  E  = 


_Q 

R2 


at  a  distance  R  from  its  center 


If  a  water  molecule  of  polarizability  a  is  placed  at  R  it 

will  be  polarized,  assuming  an  induced  dipole  moment  of 

q,q 

magnitude  y  =  aE  =  .  There  will  thus  be  an  attractive 

R~ 

force  between  the  positively  charged  ion  and  the  induced 

dipole  of  moment  y,  which  is  equal  to:  F  =  y^f  =  2y-~-  = 

2  02  dR  R 

.  The  potential  energy  describing  this  attraction  is 

related  to  the  force  F  by  Pot.  Energy  =  - / FdR  and  is  there- 

2  4 

fore  equal  to:  Pot.  Energy  =  -aQ  /2R  .  In  the  present 
calculation  of  the  above  potential  energy  which  we  call 
EPOL  is  expressed  by  equation  4.3 


ion 


n  P  (334)  (Q 
4 . 3  EPOL  =  l  l  - 

j=l  i=l  1,3 


v  2 

)  a . 


1  r  3 


2R? 


0.  is  the  ionic  charge,  a.  .is  the  polarizability  of  the 
vion  i,3 

ith  atom  of  the  jth  molecule  with  p  atoms  per  molecule  and 
n  molecules  per  cluster.  The  polarizability  of  oxygen  and 


. 
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°  3 

hydrogen  were  taken  as  0.60  and  0.44  A  respectively  (63) . 
R.  .  is  the  ion-atom  distance. 


if] 

C .  Ion-ligand  and  ligand-ligand  interaction  due  to 
dispersion  forces  (EDIS) 

The  interaction  energies  due  to  dispersion  forces 
between  ion  and  water  molecules  and  between  water  mole¬ 
cules  themselves  are  calculated  from  London's  equation  4.4 


o  u  V-,  v0 

.  ,  ^  3  h  ,  1  2  . 

4.4  E  -  -  ~6  (  )  a]_a2 

2  R  V -J-+V2 


where  v's  are  the  frequencies  of  the  ionization  limit  for 
the  ion  and  the  molecule  (61),  hv  thus  being  equal  to  the 
ionization  energy  of  the  corresponding  ion  and  the  molecule. 
The  a's  are  static  polarizabilities.  If  we  substitute  hv 
by  the  ionization  potential  in  the  above  equation,  we  obtain 
equation  4.5  used  in  the  present  calculation.  This  equation 
has  been  also  used  by  Muirhead-Gould  and  Laidler  (61)  for 
similar  evaluations. 


4.5 

a.  is  the 
1 

potential 
molecules 
cules  was 


EDIS  =  - 
2 

(point  dipole) 
of  molecule  i . 

The  distance 
taken  as  equal 


a  .  a  .  I  .  I  . 

13  13 

6 


R. 
if  3 


I  .  +1  . 
1  3 


polarizability  and  1^  the  ionization 

R.  .  is  the  distance  between  the 
1  /  3 

between  the  ion  and  the  water  mole- 


to  the  distance  from  the  ion  nucleus 


■ 
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to  the  oxygen  nucleus  (Rp  plus  0.3  A.  The  distances 
between  the  water  molecules  themselves  were  measured  from 
centre  to  centre  assuming  that  the  centre  of  the  molecules 

O 

lies  on  the  bisector  of  the  OH„  angle  and  0.3  A  from  the 

oxygen  nucleus.  The  total  dispersion  energy  was  obtained 

by  summing  up  the  energies  of  all  appropriate  molecules. 

The  ionization  potential  and  polarizability  of  water  were 

°  3 

taken  as  12.56  ev  (61)  and  1.48  A  (61)  respectively.  The 
values  of  ionization  potentials  and  polarizabilities  of 
alkali  ions  are  given  in  Table  4.1. 

D .  Dipole-dipole  repulsion  energies  (RDIP) 

The  dipole-dipole  repulsion  energy  can  be  divided  into 
three  components:  permanent  dipole  -  permanent  dipole 
repulsion  (R  ) ,  induced  dipole-induced  dipole  repulsion 
(R  )  and  induced  dipole-permanent  dipole  repulsion  (RT  ) 

roip  “  A.P.  +  RI.I.  +  RI.P. 

For  the  calculation  of  the  term  due  to  permanent  dipole- 
permanent  dipole,  R  repulsion,  the  point  charge  model 

for  the  water  molecule  was  used.  However  the  point  charge 
model  becomes  too  cumbersome  in  the  calculation  of  the 
permanent  dipole-induced  dipole  (RI  p  )  and  induced  dipole- 
induced  dipole  (R  )  interactions.  Therefore  these  were 

JL  •  -L  • 

calculated  by  using  the  point  dipole  and  point  polarizability 


of  water. 


::  I  I  •  < 
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TABLE  4 . 1 

Ionization  Potentials  and  Polarizabilities  of 

Alkali 

Ions 

Ion 

Ionization 

Potential  (eV) 

Polarizability 

°  3 
(AJ) 

Li+ 

75 . 3 

0 .03 

Na+ 

47.1 

0 .24 

K+ 

31.7 

0 . 89 

Rb+ 

27.4 

1.81 

Cs  + 

23.4 

2 . 79 

Ionization 

potentials  anc^ 

Polarizabilities,  ref  (61). 
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The  following  equations,  taken  from  references  (52)  (71) 

were  used  to  calculate  various  contributions  to  the  total 
dipole-dipole  repulsion  energy  (RDIP) 


n  1-1 


P 


4-6  Rp  p  =  I  I  l  l 

1=2  j=l  i=l  k=l 


(334)  (Q.  .)  (Q  ) 

-L  f  J  K  ^  J- 


R. 


i , j 1 


j  is  the  charge  on  the  i^1  atom  of  the  jtxl  molecule  with 

p  atoms  per  molecule  and  n  molecules  per  cluster.  R.  .  ,  , 

1 f ] ; K, I 

is  the  distance  between  the  i^*1  atom  of  molecule  and 

k^31  atom  of  the  1^  molecule. 


11  1  n  14.39  7cq  e2 

4.7  R,  1  =  I  I  1  M 

*  *  j=l  m= j+1  R 


.  S3 
ion,  3  3  ,m 


(1  +  cos  0  .  ) 

3  ,m' 


4.8  R 


n-1 


n 


I  .P  . 


1  l 

j=l  m= j+1 


2  8 . 79  4Pa..e 
M 


R2 


ion ,  3 


S3 


.m 


(1  +  cos' 


J  ,m 


a,,  is  the  polarizability  of  the  water  molecule  and  P  is  its 
M 

permanent  dipole  moment.  e  is  the  electronic  charge.  S.  is 

the  distance  of  the  molecule  from  the  m^  molecule  in 

the  cluster  and  R.  .  is  the  distance  between  the  central 

ion ,  3 

th 

ion  and  the  j  molecule.  The  numerical  factors  in  the 
above  equations  are  adjusted  to  give  the  energy  factors  in 
kcal/mole ,  for  values  of  a  in  A  ,  permanent  dipole  P  in  Debye, 
the  electronic  charge  e  in  e.s.u.  and  the  distances  R  and  S 

O 

in  A,  respectively. 


I 
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E .  Ion-ligand  molecule  repulsive  energies  (REL) 

The  ion-molecule  repulsive  forces  may  be  treated  in 

several  ways .  For  some  purposes  a  hard  sphere  model  is 

adequate,  but  it  is  somewhat  more  satisfactory  to  interpret 

repulsive  energies  by  using  a  term  of  the  form  AR  n;  A  and  n 

are  constants  and  R  is  the  internuclear  separation.  A  value 

of  12  for  n  is  assumed  in  the  well  known  Lennard-Jones 

-12 

potential.  A  term  of  the  form  AR  '  has  been  used  by  Laidler 
(6l)  in  electrostatic  potential  energy  calculations.  The 
same  form  is  adopted  in  present  calculations,  represented  by 
equation  4 . 9 

n  p 

4.9  REL  =  l  l  A.  ,R“12 
j=l  i=l 

The  evaluation  of  the  repulsion  constant  A  is  given  in  the 
next  section. 

For  the  clusters  with  more  than  one  water  molecule 
many  different  structures  are  possible.  In  the  present 
calculations  the  most  symmetrical  ligand  arrangement  was 

assumed,  where  all  m"*~ -  Ot^  (^is^:anc:es  are  equal.  In  all 

cases  the  water  molecules  were  considered  to  be  oriented 
with  respect  to  the  ion  in  such  a  way  that  the  ion  was  on 
the  bisector  of  the  H — 0 — H  angle.  The  orientation  of  the 
molecules  with  respect  to  one  another  was  determined  by 
minimizing  the  total  potential  energy  with  respect  to 
rotations  around  the  H — 0 — H  bisector. 
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The  clusters  with  two  water  molecules  were  assumed 

to  be  linear  with  the  positive  M  ion  placed  equrdistant 

between  two  water  molecules .  For  three  water  molecules , 

a  trigonal  planar  structure  was  used.  A  tetrahedron,  a 

trigonal  bipyramid  and  an  octahedron  were  used  for  four, 

+ 

five  and  six  water  molecules  bonded  to  the  ion  M  .  The 

+ 

potential  energies  E  ,  for  the  cluster  M  (H2°)n'  were 
computed  over  the  range  of  M+— — (I^O)^  distances  by  means 
of  an  IBM  360/67  computer.  The  program  used  was  taken 
from  the  work  of  S.  K.  Searles  ( 50 )  and  adapted  for  the 
present  calculations.  The  program  used  can  be  obtained 
from  the  author  on  request. 

4 . 3  Evaluation  of  the  Numerical  Value  for  the  Constant  A 

-12 

in  the  Electronic  Repulsion  Terra  REL  -  AR 
The  selection  of  the  numerical  value  for  the  constant 
A  in  the  term  REL  =  AR  '  representing  the  electronic 
repulsions  is  probably  the  most  critical  step  in  electro¬ 
static  calculations  of  the  type  used  here.  Thus,  the 
nature  of  the  central  ion  Li  or  Na  ,  etc.  xinds  expression 
in  the  calculation  essentially  only  through  the  value  of 
A.  The  dispersion  terms  also  depend  on  the  nature  of  the 
ion,  but  these  terms  are  quite  small.  For  determination 
of  the  values  of  the  constant  A  the  following  procedure  was 
adopted.  The  values  of  A  for  the  repulsions  between 
Ne...Ne,  Ar...Ar  etc.,  were  calculated  from  the  six-twelve 


. 
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potential  function  of  Lennar d- Jones  using  the  equation 
4.10,  A  =  4ea^,  e  is  the  depth  of  the  potential  well  and 
a  is  the  interatomic  distance  for  zero  potential  energy. 
The  values  for  e  and  a,  taken  from  the  literature  (64), 
as  well  as  values  of  A  calculated  from  the  equation  4.10 


are  presented  in  Table  4.2.  It  has  been  suggested  (65) 
that  the  A  values  for  the  interactions  between  two  non¬ 


identical  noble  gas  atoms  like  Ne...Ar  can  be  approximated 


quite  well  by  taking  the  geometric  mean  of  the  A  values 

p  P 

for  Ne .  .  .  Ne  and  Ar.  .  .Ar.  Thas  Ape  ^r  ~  ^  Ne,Ne  X  ^  Ar,Ar 
This  property  was  the  basis  of  the  selection  of  A  values 
for  the  present  calculation.  The  plot  of  Figure  4.2 


shows  the  values  of  A2  for  the  inert  gases  on  a  logar¬ 
ithmic  scale  versus  the  respective  Van  der  Waals  radii  of 
these  gases.  The  straight  line  obtained  allows  one  to 
select  an  A2  value  for  any  species  for  which  an  appropri¬ 
ate  radius  is  available.  It  is  important  to  realize  that 
the  selected  radius  must  be  consistent  with  the  Van  der 
Waals  radii  used  for  the  inert  gases.  Thus  the  radius 
for  the  isoelectronic  positive  alkali  ion  should  be  smaller 
than  that  of  the  corresponding  inert  gas  atom. 

One  set  of  A2  values  for  the  alkali  ion-water  repul¬ 
sions  were  obtained  by  using  the  Van  der  Waals  radii  for 
the  alkali  ions  obtained  by  R.  H.  Stokes  (66)  who  determined 


these  radii  by  quantum  mechanical  scaling  of  the  Van  der 
Waals  radii  of  the  corresponding  isoelectronic  inert  gas 


' 


. 
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TABLE  4.2 

Ip 

The  Coefficient  "A11  for  Repulsion  Term  (A/R  ~)  between 

12 

Two  Atoms  of  the  Inert  Gases  from  Equation  A  =  4ea  . 


He 1 i um-He 1 ium 

Neon-Neon 


Argon-Argon 
Krypton -Krypton 

Xenon-Xenon 


£ 

Electron-volt 
0.70  x  10”3 

3.00  x  10-3 

1.04  x  10“2 

1.41  x  10~2 

1.88  x  10~2 


a  A 

A0  kcal  -  A0 12 

2.59  0.5959  x  104 

2.77  5.5901  x  104 

3.40  3.5978  x  105 

3.60  6.1486  x  106 

4.03  3.2113  x  107 


c  is  the  depth  of  the  potential  well 

a  is  the  inter-atomic  distance  for  zero  potential  energy. 


FIGURE  4.2 


Plot  of  A2  on  logarithmic  scale  versus  the 
Van  der  Waals  radius  of  the  inert  gas  atoms. 
A  is  the  constant  in  the  repulsion  energy 
term  REL  -  AR  .  Points  shown  on  the  graph 
represent  A2  values  for  repulsion  between 
the  inert  gas  atoms:  He ---He,  Ne - Ne,  etc. 
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FIGURE  4.2 
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h 

atoms.  In  one  set  of  calculations  the  A"  term  represent- 
ing  the  water  molecules  was  taken  as  the  A2  term  of  the 

O 

oxygen  atom.  The  Van  der  Waals  radius  of  oxygen,  1.4  A, 

-12 

was  used.  The  radius  in  the  AR  term  corresponded  m 
this  case  to  ,  the  distance  between  the  alkali  and 

O  ^  i. 

"I” 

oxygen  nuclei.  The  resulting  A  values  for  the  M  ...0 
repulsions  are  given  in  Table  4.3.  The  M  . . .H  repulsions 
were  neglected. 

A  second  set  of  A  values  was  obtained  with  the  A 

h 

values  of  the  alkali  ions  as  determined  above  and  the  A 

value  for  water-water  molecule  repulsions  given  in  the 

-12 

literature  (p.214  ref  (64)).  The  radius  in  the  AR  term 

o 

in  this  case  corresponded  to  R^  -^  +  0.3  A. 

A  third  set  of  A  values  was  evaluated  for  Na  (OHpJn 
in  which  the  A  value  was  taken  to  be  the  same  as  that  for 
the  Ne . . . Ne  repulsions.  The  justification  for  this  pro¬ 
cedure  lies  in  the  observation  that  Na+  should  be  somewhat 
smaller  than  Ne  while  H20  should  be  somewhat  larger 
therefore  an  approximate  cancellation  of  the  two  errors 
might  be  expected. 

The  potential  energies ,  calculated  by  using  the 
first  set  of  A  values  are  given  in  Table  4.3  for  all 

-j- 

alkali  ions.  The  potential  energies  for  Na  (H20)  calcu¬ 
lated  from  three  sets  of  A  values  are  plotted  in  Figure  4.3. 

The  above  three  choices  of  A  values  are  consistent 
with  the  accepted  set  of  Lennard— Jones  parameters  for 
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TABLE  4 . 3 


+ 


Potential  Energies  of  Clusters  M  (H20)  from  Electrostatic 


Calculations  Using  Repulsion  Term  AR 


-12 


V. 

0 

W 

<1 

1 

n 

Li+ 

Na^ 

K+ 

Rb+ 

Cs  + 

0,1 

53.3 

42.9 

31.2 

28 . 1 

23 . 1 

1 — 1 

CO 

Ph 

1.56 

1.74 

2 .04 

2 . 16 

2 . 38 

0 , 2 

95.9 

79.1 

58 . 8 

53.6 

44.4 

R3, 1 

1.60 

1.76 

2.06 

2 . 18 

2 .38 

0,3 

119 . 6 

102.7 

79 . 8 

73 . 7 

62.0 

R3 , 1 

1.66 

1.80 

2.08 

2.20 

2 . 42 

0,4 

129 . 7 

116.5 

94 . 9 

89 . 1 

76 . 3 

i — 1 

CO 

Pi 

1.72 

1.86 

2.12 

2.24 

2 .44 

0,5 

123.8 

117.4 

101.6 

97.5 

85.4 

R3 , 1 

1. 84 

1.94 

2 .18 

2 .28 

2 . 48 

0 , 6 

118.4 

116.0 

105.9 

103.6 

92 . 6 

R3,l 

1.96 

2 .04 

2.24 

2.34 

2.52 

A  x 

kcal 

105 
x  U2 

0.0298  0.0861 

0 .4200 

0 . 8400 

2 . 100 

(a) 

Repulsion  term 

determined 

according  to 

procedure 

described 

in  this 

Chapter 

from  noble  gas  A  values,  Van  der 

W  a  a  1  s 

radii  of  noble 

gases  and 

Stokes  radii 

of  M  . 

AE  in  kcal/mole 
o,  n 

R  in  Angstroms  corresponds  to  distance  from  centre  of 

3 , 1 

ion  to  nucleus  of  oxygen  atom  of  water  molecule. 
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the  inert  gases.  These  parameters  were  determined  from 
second  virial  coefficients  which  are  due  to  weak  inter- 
molecular  interactions.  Assuming  that  these  parameters 
are  correct  in  the  weak  interaction  range  they  will  be 
applicable  to  stronger  interactions  (smaller  intermolecu- 
lar  distances)  only  if  the  Lennard-- Jones  potential 
function  is  of  the  correct  form.  It  would  be  desirable 
to  obviate  the  use  of  the  Lennard- Jones  six,  twelve 
potential  by  dropping  both  these  terms  (EDIS  and  EREP) 
from  the  potential  energy  expression  and  replacing  them 
with  experimentally  determined  potentials  based  on  molecu¬ 
lar  beam  scattering  experiments.  The  scattering  in 
relatively  high  kinetic  energy  beams  is  due  to  strong 
interactions  occurring  at  close  distances  between  the 
colliding  partners.  Unfortunately,  experimental  potential 
functions  for  the  present  pairs  Li^.-.H^O,  Na*...!^!)  are 
not  available.  However,  as  shown  previously  the  inter¬ 
action  for  the  pair  Ne. . .Ne  appears  to  be  a  good  approxim¬ 
ation  for  that  between  Na+  and  II^O.  Therefore,  calculations 
were  made  for  sodium  using  the  potential  term  7190/9.99 
kcal/mole  instead  of  EDIS  and  REL.  This  potential  energy 
term  was  determined  by  Amdur  from  experimental  studies  of 
Ne...Ne  atom  collisions  at  high  energies  (67). 

The  numerical  results  calculated  on  the  basis  of  the 
assumptions  discussed  above  will  be  considered  in  the  next 


section . 


. 
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4.4.  Discussion  of  Electrostatic  Calculations 

The  potential  energies  E  ,  where  n  represents  the 

number  of  water  molecules  (n=l  to  6)  bonded  to  the  central 

ion,  were  calculated  for  all  the  alkali  ions  using  an  IBM 

360/67  computer.  The  potential  energy  minima  found  are 

given  in  Table  4.3,  where  the  notation  AE  is  used  in 

o ,  n 

analogy  with  the  notation  for  the  experimental  AHq  ^ 
values.  The  electrostatic  potential  energies  AE^  ^  cor¬ 
respond  best  to  the  heats  of  hydration  AH  at  0°K, 
because  of  the  neglect  of  all  thermal  motions  in  the 
calculations.  The  comparison  at  298°K  assumes  the  AC  =  0. 
This  assumption  does  not  involve  considerable  approxi¬ 
mation,  since  AC^  for  reaction  0,1  for  example  can  be 

P 

5 

roughly  estimated  to  be  R.  The  value  of  AC  is  estimated 

2  p 


0,1 


+  + 

M  +  H20  +  M  (H20) 


by  taking  into  account  only  the  translational  and  rotational 

motions  of  molecules.  It  was  also  assumed  that  H20  and 

+ 

M  (II20)  have  the  same  rotational  contributions  to  the 

5 

heat-capacity.  Thus  the  value  of  AC^  =  R  would  be  mainly 
due  to  the  loss  of  translational  degrees  of  freedom  by 
M+ .  Vibrational  contributions  to  the  heat-capacities  of 
H20  and  M+(H20)  were  not  taken  into  account  since  they 
make  a  very  small  contribution.  Vibrational  motion  can 
contribute  any  amount  (between  O  and  R)  to  the  heat 


. 
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capacity,  the  later  value  is  approached  only  at  high  tem¬ 
peratures  for  most  molecules  (68) ,  However,  since  the 
product  M+ ( H ^ 0 )  in  reaction  0.1  has  more  vibrational 

degrees  of  freedom  than  the  reactant  h^O,  the  absolute 

5 

value  of  ACp  will  probably  be  even  smaller  than  -jR  or 
-5  cal/degree  mole.  At  298°K  AC^AT  is  approximately  equal 
to  ~-1.5  kcal/mole 

-1 

Thus  if  we  assume  ACpAT  ~  -1.5  kcal/mole  ,  then 
by  use  of  equation  4.11  which  expresses  the  temperature 
dependence  of  reaction  heats, 


4.11 


2  9  8  °  K 

+  f  AC  dt 

J  P 

0°K 


we  obtain  the  following  relation: 


AH0,1  s  AE0,1  “  1-5  kcal 


The  AH  is  negative  value,  therefore  its  absolute  value 

0,1  -1 

should  be  larger  by  approximately  1.5  kcal/mole  than  the 

I AE^  -i  I  .  The  difference  of  1.5  kcal  between  AHn  and 
IQ,!1  u  '  x 

AE,  is  not  significant.  Therefore,  it  is  felt  that 

U  /  -L 

meaningful  conclusions  may  be  drawn  by  a  comparison  of  the 
electrostatic  calculations  with  experimental  enthalpy 
values,  without  taking  into  consideration  the  above  cor¬ 


rection  . 


Given  in  Table  4.4  are  the  values  of  the  indiv 

+  .  • 

idual  energy  terms:  EDIP ,  EPOL  etc.,  for  the  Na  ion  m 
order  to  illustrate  the  relative  magnitude  of  each  contri- 
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TABLE  4 . 4 


Individual  Energy  Contributions 

Energy  Calculations  of  Cluster 


to  Potential 

Na+  (0Ho) 

2  n 


n 

-EDIP 

-EPOL 

-EDIS 

RDIP 

REL 

1 

36 . 871 

14.982 

2 .283 

0 

11.180 

2 

72 . 243 

28.741 

4.425 

6 . 810 

19 . 494 

3 

104.064 

39 . 720 

6 .531 

25.236 

22 . 330 

4 

130 . 794 

47.015 

8.353 

49 . 490 

20 . 117 

5 

151.423 

50 . 385 

10 .680 

80 .234 

15.149 

6 

165 . 773 

50 . 341 

11.608 

104.830 

9 .946 

(a)  All  values  in  kcal/mole. 
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bution.  Figures  4.3  and  4.4  give  plots  of  AE  and  the 

O  f  1 1 

experimentally  measured  AHq  n  for  sodium  and  cesium.  It 
can  be  seen  from  the  Figures  that  for  all  these  different 
choices  of  "A"  value,  the  calculated  values  are  consider¬ 
ably  larger  than  the  experimental  values.  In  general,  the 
calculated  AEQ  ±  and  AEQ  2  are  about  70%  larger  than  the 
corresponding  enthalpy  values. 

At  present  we  cannot  understand  very  well  why  all 
choices  of  A  values  lead  to  AE  results  which  are  too 

\J  /if 

high.  See  Figure  4.3.  Figure  4.3  also  shows ,  that  AEQ,n 
values,  calculated  from  the  Amdur  empirical  term  719Q/r9.99 
kcal/mole,  are  approximately  as  high  as  the  AEQ,n, 
calculated  from  different  choices  of  A  values  derived 
from  the  Lennard-Jones  potential  function.  Thus  the 
Amdur  empirical  term  also  represents  a  "soft"  repulsion 
between  the  ion  and  molecule.  The  following  consideration 
may  be  an  explanation  for  the  differences  between  calcu¬ 
lated  and  experimental  values .  All  choices  of  A  depend 
on  a  certain  assumed  size  parameter  (radius)  for  the 
positive  ion  and  the  water  molecule.  However,  since  the 
electron  cloud  of  the  positive  ion  is  compressed,  because  of 
PHe  larger  nuclear  charge,  the  electronic  repulsive  forces  in 
a  strong  collision  between  the  positive  ion  and  the  neutral 
should  be  larger  than  those  occurring  in  a  collision  of  two 
neutrals  of  the  same  size.  This  stiffness  of  the  electron 
cloud  of  the  alkali  ion  is  not  taken  into  account  in  the 
present  choices  of  A.  An  empirical  way  of  correcting  for 


FIGURE  4.3 


Comparison  between  experimental  enthalpies  AHQ,n 

and  potential  energies  AEq  n  obtained  on  the 

basis  of  electrostatic  type  calculations  for  sodium 

hydrates  Na+(H20)n.  0  experimental  AHQ,n,  ©  AE0,n 

calculated  by  adjusting  value  of  constant  A  in 

repulsion  term  AR  so  as  to  make  AH^  =  AE^  ^ 

01  AE  with  A  value  for  Na+  —  "0  repulsion 
o  ,n 

+ 

obtained  from  Stokes  radius  for  Na  and  Van  der  Waals 
radius  of  oxygen,  □  calculated  with  A  value  for 

Na+ - H20  repulsion  set  equal  to  A  value  for 

Ne-Ne  repulsion  (65)  A  calculated  with  A  value  for 

"4“  .  -t" 

Na  - - H20  repulsion  based  on  Stokes  radius  of  Na 

and  A  value  for  H-^O-IUO  repulsion  (64),  ®  calculated 

with  the  potential  term  7190/^9.99  kcal/mole 

instead  of  EDIS  and  REL  (67)  0  E  calculated 

o  f  n 

by  Garrick  (71) . 
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FIGURE  4.4 


Comparison  between  experimental  enthalpies  AHQ^n 

and  potential  energies  AEq  n  obtained  on  the 

basis  of  electrostatic  type  calculations  for 

cesium  hydrates  Cs+(H20)n.  0  experimental 

AH  ,  ©  AE  calculated  by  adjusting  the 

o ,  n  o,n 

-12 

value  of  constant  A  in  repulsion  term  AR 
so  as  to  make  AH^  ^  =  AEQ  ©  AEQ^n  with  A 

value  for  Cs+- — 0  repulsion  obtained  from 
Stokes  radius  for  Cs+  and  Van  der  Waals  radius 
of  oxygen,  0  calculated  with  A  value  for 
Cs+-H20  repulsion  based  on  Stokes  radius  of 
Cs+  and  A  value  for  H20-H20  repulsion,  A 
calculated  with  A  value  for  Cs+-H20  repulsion 

set  equal  to  A  value  for  Xe-Ne  repulsion, 

E  calculated  by  Garrick  (71) . 

U  o,n 
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the  greater  "stiffness"  of  the  positive  ion  electron  cloud 
is  to  choose  a  bigger  radius.  For  example  the  radii  of 
Pauling  (69)  or  Gourary  and  Adrian  (70)  might  have  been 
selected.  The  difficulty  in  such  a  procedure  lies  in 
the  correlation  of  the  ionic  radii  with  the  radii  of  the 
inert  gases,  which  are  the  bases  of  the  plot  in  Figure 
4.2.  This  becomes  more  obvious  when  the  Na  radius  obtained 
by  Gourary  and  Adrian  is  compared  with  the  atomic  radius  of 

neon  and  it  is  found  that  the  radius  of  Na+  is  bigger 

,  ° 
than  that  of  Ne  (Na  =  1.17,  r  -  =  1.16  and  Ne  =  1.12  A). 

This  is  not  of  course  acceptable,  and  shows  that 
the  Gourary  and  Adrian  ionic  radii  are  probably  too  low 
for  negative  ions  and  too  high  for  positive  ions  by 
approximately  10%  (3Q  )  • 

4 . 5  Comparison  of  Electrostatic  Potential  Energies  with 

Previous  Calculations . 

The  present  calculations  are  very  similar  to  those 
of  Garrick  (71)  who  used  a  point  dipole  model  to  cal¬ 
culate  the  potential  energies  of  a  number  of  hydrated  ions 
including  Na+(H20)^,  Na  (H20)  K  (H^O)^,  Rb  (I^O)  ^  and 
Cs+(H20)6  respectively. 

-12 

For  the  ion-molecule  repulsion  (REL  =  AR  ) ,  the 
equation  used  by  Garrick,  (REL  =  XR_V)  was  adapted  from 
Lennard-Jones  calculations  of  inert  gas  ion-atom  repul¬ 
sions.  The  X  in  Garrick's  expression  for  REL  corresponds 
to  A  value  in  the  present  calculations.  The  values  of 


111. 


force  constant  A  and  force  index  v  for  different  ions 
were  deduced  by  means  of  certain  more  or  less  plausible 
assumptions  (71) . 

Garrick's  a  priori  calculations  for  M  (f^O)  g  gave 
the  equilibrium  distance  between  the  point  ion  and  point 
molecule  as  2.22  (2.04),  2.58  (2.24)  ,  2.71  (2.34)  ,  2.92 
(2.52)  A,  for  M+  =  Na,  K,  Rb ,  Cs ,  respectively.  In  com¬ 
parison,  the  values  for  the  ion-oxygen  atom  (R3#1)  equil¬ 
ibrium  distance,  from  this  work  are  given  in  the  brackets. 
The  values  for  EDIP  +  EPOL  +  RDIP  are  135  (115)  ,  102 
(108)  ,  93  (101)  ,  81  (92)  kcal/mole  and  REL  are  18  (10)  , 

18  (16) ,  15  (18) ,  14  (19)  kcal/mole  for  hexahydrate  of 

Na,  K,  Rb  and  Cs  respectively.  The  total  potential 
energies  in  Garrick's  calculations  are  114  (116),  83  (105), 

78  (103) ,  67  (92)  kcal/mole,  compared  to  present  calcu¬ 

lations  in  the  brackets.  The  higher  total  potential 
energies  obtained  here  are  mainly  due  to  dispersion  forces 
^j^pch  were  not  included  in  Garrick  s  calculations . 

Besides  Garrick/  Laidler  (72)  has  performed  point 
dipole  calculations  for  a  number  of  hydrated  alkali  ions. 
Unfortunately  most  of  the  M+(H20)r  results  were  presented 
as  the  total  of  a  cycle  dealing  with  the  heat  of  hydration 
of  the  ion.  The  electrostatic  energy  of  formation  for 
Li+(H20)^  and  Rb+(H20)^  obtained  by  Laidler  were  120  and  67 
kcal/mole,  compared  to  129  and  89  kcal/mole  in  present 


calculations . 
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Both  Laidler  and  this  author  used  the  same  depend- 

-12 

ence  of  the  ion-molecule  repulsion  (REL  =  AR  )  ,  but 

Laidler  adopted  a  semi-empirical  method  to  evaluate  the 

constant  A  in  a  repulsion  term.  The  value  of  A  was 

adjusted  so  as  to  give  a  distance  between  the  ion  and  the 

water  molecule  equal  to  the  sum  of  the  Goldsmith  ionic 

radius  and  the  radius  for  the  water  molecule.  Laidler* s 

value  for  Rb+(H20)4  is  closer  to  our  experimental  value 

AH  =53.0  kcal/mole ,  than  to  our  calculated  one.  The 
0/4 

values  of  the  constant  A  and  the  equilibrium  distances 
between  the  ion  and  the  water  molecules  for  Rb  (1^0)4  or 
Li+(Ho0)A  are  not  reported  in  Laidler' s  work,  so  that  a 
more  detailed  quantitative  comparison  of  the  two  calcu¬ 
lations  cannot  be  made. 

O 

4.6  The  Comparison  of  Calculated  ^En_lfTl  with  Experi¬ 
mental  for  M+(H20)n_I  +  H20  ->  M+(H2Qn 

One  of  the  purposes  of  doing  electrostatic  calcu¬ 
lations  was  to  make  a  qualitative  comparison  between 
predicted  and  experimentally  obtained  enthalpy  changes  for 
M+(H20)n  complexes.  Since  the  calculated  AEQ,n  obtained 
were  found  to  be  quite  different  from  the  experimental 
results  we  decided  to  utilize  the  calculations  only  for  a 

O 

comparison  of  the  relative  changes  of  AEn_^n  predicted 

O 

by  calculation  and  experimentally  determined 

calculate  such  relative  changes  it  appeared  most  meaningful 


. 
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to  determine  an  empirical  value  for  the  constant  A  in  the 
term  AR  ^  by  setting  ^  =  Hq  ^  .  The  evaluation  of 

value  A  by  this  method  is  given  in  the  following  paragraph. 

To  calculate  the  m'*" - repulsion  energy  from  the 

term  AR  it  is  necessary  to  assign  values  to  AQ ,  a  con¬ 
stant  for  the  ion-oxygen  atom  repulsion  and  AR ,  a  constant 
for  the  ion-hydrogen  atom  repulsion.  Minimization  of  the 
total  potential  energy  equations  for  M  (E^O)  with  respect 
to  j  leads  to  equation  4.12  from  which  AQ  can  be  cal¬ 
culated  as  a  function  or  R^  .  assuming  each  value  of  R± ^ 


corresponds  to  the  equilibrium  ion— ligand  distance.  The 
terms  containing  A^  may  be  expected  to  be  negligible  com¬ 
pared  with  the  AQ  terms.  The  Aq  is  expected  lo  be  larger 
than  Ah  since  the  repulsion  depends  on  the  number  of 
electrons  associated  with  the  atom  0J) .  Also,  the  ion- 
hydrogen  distances  are  much  greater  than  the  ion— oxygen 
distance.  Neglecting  the  A^ ,  one  obtains  a  relation  of 

the  type : 

4.13  Aq  =  f(R3fl) 

which  gives  AQ  as  a  function  of  ion-oxygen  distance,  R-^-^- 
The  potential  energy  of  the  cluster  M  (H^O)  ,  which  depends 
on  the  parameters  R^  ^  an<^  ^iven  k Y  equation  4.14 


. 
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4.14 


AEq , ±  F r3  ,  . 


In  order  to  obtain  a  value  for  AQ  and  ^  at  equilibrium 
we  equate Ae^  ^  with  the  experimental  AHq  ^  which  leads 
to  4.15 


4.15 


AH0,1  -  F(V  R3,l> 


Expressions  4.13  and  4.15  represent  a  system  of  two 
equations  with  two  unknowns  (Aq,  R3  ^  e<l*)  which  can  be 
solved.  The  values  of  R3  1  with  the  corresponding  numbers 
for  Aq  (obtained  by  4.13)  are  substituted  into  equation 
4.14.  This  procedure  leads  to  the  AEQ  ^  as  a  function  of 


R 


3  1.  The  value  of  R3  1  corresponding  to  the  AEQ  1  -  AH0  ^ 
is  R_  equilibrium  distance  between  the  centre  of  the  ion 

O  f  1. 

M+  and  the  centre  of  oxygen  atom  in  the  water  molecule. 

The  associated  value  for  A^  which  corresponds  to  R3  ^ 

equilibrium  distance  was  then  used  to  calculate  REL  and 

Ae  for  higher  clusters, 
o ,  n 

The  AE  values  obtained  by  this  procedure  were 
then  used  to  evaluate  AEn_^n  =  AEQ^n  -  AEQ^n_1  for  all 
ions,  for  the  range  of  n=l  to  6.  See  Table  4.5.  In 
Figure  4.5,  AEQ  ^  is  shown  as  a  function  of  R3  ,  1  for 
M+  +  H90  M+(H90)  .  The  energy  minimum  for  each 

individual  ion  corresponds  to  the  total  calculated  binding 
energy  of  the  clusters  relative  to  M  and  H^O  molecules 
at  infinite  distance. 
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TABLE  4.5 


Potential 

Energies 

of  Clusters 

M+ ( H  0  0 ) 

2  'n 

from  Electrostatic 

Calculations 

Using  Repulsion  Term 

AR  12  a 

-AE 

o. 

n 

r  •  ^ 

Li 

Na+ 

K+ 

Rb+ 

Cs  + 

0,1 

34 . 0 

24.0 

17.9 

15.9 

13 . 7 

i — 1 

<T) 

1.92 

2 .26 

2 .60 

2 . 76 

2.98 

0 , 2 

63.5 

45 . 5 

34.4 

30.9 

26.5 

R3, 1 

1.92 

2 .26 

2 .60 

2 . 76 

2 .98 

0,3 

85.8 

63.0 

48.7 

43 . 9 

37 . 9 

R3 , 1 

1.96 

2.30 

2.62 

2.76 

2.98 

0,4 

100 . 0 

76.7 

60 . 6 

55.0 

47.9 

R3 , 1 

2 . 02 

2.34 

2 . 66 

2 . 76 

2.98 

0,5 

105.3 

85 . 1 

69.2 

63.2 

55.9 

i — 1 

ro 

P4 

2 . 10 

2.38 

2 .68 

2 . 78 

3.00 

0,6 

107.2 

90  .  8 

76 . 0 

70.5 

63.5 

i — 1 

00 

Ph 

2.16 

2 .44 

2 . 74 

2.84 

3.06 

A  x 

5 

10 

°  12 
x  A 

0 . 19 

0.95 

3.9 

7.37 

15 . 87 

kcal 

(a)  A  values  for  term  AR^  ^  in  electrostatic  calculations 
obtained  by  setting  AH^  ^  =  AE^  ^ 

O 

R_  in  A  corresponds  to  distance  from  centre  of  ion  to 

^  r  -L 

nucleus  of  oxygen  atom  of  water  molecule. 


Potential  Energy  E(  kcal  / mole) 
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R  3,1  (^) 

Plot  of  the  potential  energies  evolved  for 
the  reaction  m"*~  +  H^O  ->  M  (H2O) 


FIGURE  4.5 


FIGURE  4.6 

_ _ — - —  O 

Comparison  of  experimental  AHn_1  n  and 

O 

potential  energy  differences  AEn_^n, 
ion  clusters  0  AH  of  Li  , 

□  Ah  of  Na+,E1  AE  of  Na+ ,  A  AH  of  Cs  + 


calculated 
for  alkali 
©  AE  of  Li+, 
,  A  AE  of 
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FIGURE  4.6 
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o  O 

Figure  4.6  shows  a  plot  of  AHn„j_  n  and  AEn-l,n  aS  a 
function  of  n,  for  the  small  ion  Li  ,  the  medium  sized  ion 
Na+  and  the  large  ion  Cs+ .  According  to  the  theoretical 

O 

model,  the  decrease  of  AE  with  n  is  due  to  an 

increase  in  dipole— dipole  repulsions .  As  can  be  seen  from 

Figure  4.6  the  electrostatic  calculations  result  in  values 

of  AE °  which  qualitatively  agree  with  the  experimental 

n  jl  f  n 

o  ° 

AH  ,  .  Thus  the  changes  in  AH  ,  should  be  mostly 

n  - 1 ,  n  n  -  ±  ,  n 

due  to  increases  in  dipole-dipole  repulsions.  However, 
for  small  n  (n=2  and  3)  it  is  found  that  AEr_-^r  falls  off 
less  rapidly  than  the  AH^-j^  .  This  effect  is  most  pro¬ 
nounced  for  Li+  and  is  almost  unnoticeable  for  Cs  .  ihe 


more  rapid  initial  drop  of  the  AHr_-^r  terms  could  be 
due  to  the  presence  of  chemical  (non-electrostatic)  bonding 
in  the  M+OH  complex.  This  bonding  should  occur  only 


for  the  attachment  of  the  first  water  molecule  and  be  com 


pletely  or  almost  completely  absent  for  the  other  water 
molecules.  Such  a  state  of  affairs  would  exist  if  a 
certain  amount  of  dative  bonding  between  the  oxygen  lone 
pair  and  the  lowest  empty  orbital  of  the  alkali  ion 
occurred.  The  oxygen  lone  pair  electrons  would  then 
partly  screen  the  positive  charge  from  the  other  water 
molecule  (73) .  The  effect  should  be  largest  in  lithium 
which  has  the  most  covalent  character  and  should  be  weakest 
with  cesium.  This  is  what  is  observed  in  Figure  4.6. 

At  higher  n  the  AE^  R  curve  is  found  to  decrease  more 
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rapidly  and  cross  the  AH  ,  ,  n  curve.  For  Li  and  Na  , 

n.  _L  f  n. 

AE .  q  and  AE  are  very  much  smaller  than  the  corresponding 
AH  (Figure  4.6) .  As  mentioned  before  the  calculations  are 
made  under  the  assumption  that  all  water  molecules  are 
equidistant  from  the  ion.  For  Li  and  Na  which  are  small 
ions  this  assumption  is  probably  inadequate  at  n=5  and  n=6. 
Thus  the  cross-over  at  n=4  should  indicate  that  at  this 
point  a  change  to  a  new  structure  occurs  in  which  all  water 
molecules  are  not  equivalent.  The  AE^  ^  and  AE^g  for  Cs 
are  somewhat  smaller  than  the  corresponding  AH  values . 
However,  the  differences  between  the  calculated  and  experi¬ 
mental  values  are  very  small  so  that  one  may  conclude 
that  the  data  do  not  indicate  a  change  of  structure  at 
n=4  or  5  but  a  persistence  of  the  equivalence  of  all  water 
molecules  to  about  n=6 .  This  is,  of  course,  consistent 

with  the  larger  radius  of  Cs+ 

Yu.  A.  Kruglyak  (74)  investigated  the  theory  of 
the  solvation  of  ions  having  external  completed  sp-shells. 
He  suggests  that  the  nature  of  bonds  in  ion-solvates 
interactions  can  be  understood  only  on  a  quantum  mechani¬ 
cal  basis;  purely  electrostatic  interpretations  are  not 
sufficient  and  many  facts  can  only  be  explained  by  means 
of  donor-acceptor  bond  representations.  This  supports  our 
evidence  of  the  presence  of  chemical  (non-electrostatic) 
bonding  in  the  M+OH2  complex,  explained  previously. 


' 
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4.7  Calculation  of  the  Entropy  of  M1 (H^O) 

The  experimental  entropy  change  for  a  reaction 
such  as  M4  +  H^O  ->  M+(H20)  can  be  compared  to  the  theor¬ 
etical  value  obtained  by  means  of  statistical  mechanics. 

The  change  of  the  entropy  for  such  a  reaction  can  be  repre¬ 
sented  as  the  sum  of  translational,  rotational  and  vib¬ 
rational  entropy  charges. 


4.14  AS  =  AS t  +  ASr  +  ASv 


Each  term  of  the  right  hand  side  of  equation  4.14  is  then 


evaluated  separately. 

The  translational  entropy  change  for  the  condensa¬ 


tion  reactions  considered  here  is  given  by  the  Sackur' 


Tetrode  equation  4.15  (75) 
4.15 


3  NM+ (H20) 

ASt  =  2R  Ln  N.._l  x  N 


R  In  T  +  2.311 


h2o 


where  Njyi+(H  O)  '  an(^  NH  0  are  mo-Lecu^-ar  weights  (g/mole) 

2  2 

of  m+,  ,  M+  and  H00  respectively.  T  is  the  absolute 

(H20)  2 

temperature  in  degrees  K;  the  gas  constant  R  must  be 


substituted  in  cal/mole— degree  to  obtain  AS^_  in  cai/mole 
deg . 

The  rotational  entropy  change  can  be  calculated 
from  the  rotational  entropy  of  H20  and  (H20) .  The 
appropriate  equation  (75)  is  4.16 

4.16  S°  =  R(ln  Qrot  +  D 
rot 
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where  Q  ,  the  rotational  partition  function  (75)  is  given 
by  4.17 


4.17 


Q 


(8)  5//2  (KT)  3/2  (III  )  1/7  2 

x  y  z 


rot 


ah 


3 


where  I  ,  I  ,  I  are  the  three  moments  of  inertia  about  the 

x'  y  z 

three  principal  axes  and  "a"  is  the  symmetry  number.  When 

equation  4.16  is  applied  to  calculate  the  change  of  entropy 

+  + 

for  the  reaction  of  M  +  H20  ^  (H20)  r  the  following 

equation  will  be  obtained 


4.18 


o  [Qrot]M+(H  0) 

s  =  r  In  — — - - — 

r0t  [Q  .]H90 

L^rot  2 


When  the  proper  values  of  the  partition  functions  are 
substituted  in  equation  4.18  and  appropriate  simplification 


has  been  made,  the  S  ,  becomes 

rot 


,+ 


4.19 


_  R  tWz]M  (H20) 

>rot  2 

[III  ] Ho0 
x  y  z  2 


The  moments  of  inertia  of  H20  are  (76)  Ix  =  1.92,  I  =  1.02 

and  I  =  2.95  x  10_4°  g  cm2.  The  moments  of  inertia  for 
z 

1  + 

the  cluster  M  (H20)  are  based  on  the  M  . . .0H2  equilibrium 
distances.  As  shown  in  section  4.7  reasonable  distances 
between  the  centers  of  the  M+  ion  and  the  oxygen  atoms 

were  obtained  from  electrostatic  calculations.  The  M - 0I12 

equilibrium  distances  as  well  as  moments  of  inertia  are  given 
in  Table  4.6.  Substitution  of  the  moment  of  inertia  into  4.19 
leads  to  rotational  entropy  change,  given  in  Table  4.9. 
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TABLE  4 . 6 


Moments  of 

Inertia  and  M+-OH2 

Equilibrium  Distances 

Cluster 

m+-oh2 

Distance 

O 

A 

I  X 

4 

X 

in40 

10  g  - 

I 

Y 

2 

cm 

z 

Li+ (H20) 

1.92 

1.95 

33.2 

35.0 

Na+ (H00 

2.26 

1.95 

87.3 

89 . 0 

K+ (H20) 

2 . 60 

1.95 

147.0 

149.0 

Rb+  (I4?0) 

2.76 

1.95 

182 . 4 

184 . 0 

Cs+  (H20) 

2.98 

1.95 

233.0 

235.0 

■ 
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The  z  axis  is  mutually  perpendicular  to  the  x 
and  passes  through  the  point  (0.04,  0,  0) 


Principal  axes  of  Li+ (H^O) 

The  z  axis  is  mutually  perpendicular  to  the  x 
and  passes  through  the  point  (0.02,  0,0) 


FIGURE  4 . 7  Principal  axes  of  I^O  and  Li  (H^ 


0) 

and  y  axes 

.22,  0.77,  0) 

X 

.22,  -0.77,  0) 

and  y  axes 

O)  used  to 


calculate  the  moments  of  inertia. 
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The  three  principal  axes  for  the  H20  and  Li+(H20)  used  to 
calculate  the  moments  of  inertia  are  shown  in  Figure  4.7. 

To  calculate  the  vibrational  contribution  to  the 
entropy  change,  it  was  assumed  that  the  fundamental  vib¬ 
rations  of  water  molecule  remain  the  same  before  and  after 
the  formation  of  the  cluster  M+(H20).  The  coordinated 
H20  is  treated  as  a  rigid  molecule  which  has  a  stretching 
vibration  along  the  M  — — —oxygen  axis  and  two  rocking 
vibrations  about  its  y  and  z  axes. 

In  the  calculation  of  the  stretching  vibrational 
frequency,  the  water  molecule  and  the  M  ion  were  treated 
as  two  mass  points.  The  force  constant  for  this  vibration 
was  estimated  by  fitting  an  equation  of  a  parabola  to  the 
bottom  portion  of  the  potential  energy  curves  of  Figure  4.5 

obtained  in  section  4.7.  The  frequency  v  then  can  be  cal- 

1  k  t 

culated  from  the  equation  4.20  v  =  7=; — (— )  2  where  u  is  the 

^  2  tt  y 

reduced  mass  of  v/ater  — M*1  ion  system  and  k  is  the  force 

O 

constant.  From  e  =  hv  and  equation  4.21  (75)  s  . .  were 

o  a  vrb 

calculated.  In  Table  4.7  the  force  constants  k,  the  vib¬ 


rational  frequencies  v  and  entropy  contribution 
are  tabulated. 


S 


stretching 


O 


s 


vib 


R  e 

( - )  (—  ) -R  In  [1-exp  (-e  /RT)  ] 

-1  Ft  0 

exp (e  /kT) 


4.21 


. 
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TABLE  4 . 7 

Force  Constants,  Vibrational  Frequencies  and  Entropy 
Contributions  Due  to  Stretching  Vibrations  of 


m+-oh2 

Cluster 

Force  Constant 

-4 

k  x  10  dyne/cm 

Frequency 

n  12  -1 
v  x  10  sec 

Entropy 

AS  ° 

a  ( 2 9  8  °K) e . u . 

Li+ (H20) 

11.70 

18.0 

0.50 

Na+ (H20) 

8.10 

11.05 

0 .60 

K+ (H20) 

5 . 80 

8.47 

1.60 

Rb+  (HO) 

4 . 35 

6 .68 

2.00 

Cs+  (H20) 

3.50 

5 . 80 

2.27 
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The  two  rocking  vibrational  frequencies  can  be 
approximately  evaluated  from  equation  4.22  (76)  where  n  is 
the  symmetry  number,  AP.E.  is  the  potential  energy  barrier 
for  rotation,  and  I  is  the  moment  of  inertia. 


2  21 

AP.E.  can  be  calculated  from  the  coordinates  of  the  molecule 

in  its  minimum  energy  position  and  its  maximum 

energy  position  (rotation  by  180°  about  the  y  and  z  axes 

from  the  minimuni  position)  by  use  of  equations  4.2,  4.3,  4.5 

O 

and  4.9.  From  e q(  }  =  hv(y,z)  and  ec2uation  4’21'  Sy  and 

S  were  calculated.  In  Table  4.8  AP.E  and  AP.t^,,  the 
z 

vibrational  frequencies  and  v z  and  entropy  contribution 

o  o 

S  and  S  are  tabulated, 

y  z 

The  total  entropy  change  can  now  be  calculated 
by  summing  the  individual  entropy  contributions .  As  can 
be  seen  from  Table  4.9  relatively  good  agreement  is  obtained 
between  the  calculated  and  observed  entropy  changes. 

Entropy  values  calculated  at  298°K  and  entropy  values 
calculated  at  the  temperatures  corresponding  to  the  middle 
of  the  van't  Hoff  plot  temperature  range  from  which  the 
experimental  AS  values  were  obtained  are  shown  in  Table  4.9. 
The  entropy  changes  for  the  formation  of  the  larger 

clusters  [M+(H20)n  for  n  9reater  than  one^  can  also  be  cal¬ 
culated  but  the  above  procedure  would  involve  an  increasing 
number  of  assumptions.  Therefore,  such  calculations  were 
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not  performed. 

It  is  interesting  to  note  that  the  experimental 

-AS  increases  with  increase  of  n.  The  same  result  is 

n-1,  n 

obtained  for  the  halide  ions  by  M.  R.  Arshadi,  Ph.D.  thesis. 
The  increase  of  the  AS  values  occurs  even  though  the  bonding 
forces  (represented  for  instance  by  AHn_^n)  decrease.  It 
is  reasonable  to  assume  that  the  increase  of  the  AS  terms 
reflect  the  gradual  loss  of  freedom  due  to  crowding  of 
the  molecules  around  the  central  ion.  The  study  of  the 

O 

H+(H  0)  clusters  (24)  showed  that  the  ASn_lin  values 
decrease  somewhat  with  increase  of  n.  This  difference  of 
behaviour  between  the  clusters  around  the  centrally  located 
(large)  ions  and  water  clusters  containing  an  extra  proton, 
may  be  considered  as  evidence  for  the  mobility  or  non¬ 
localization  of  the  charge  in  the  hydrates  H  (H20)n* 
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5. 


APPLICATION  OF  EXPERIMENTAL  RESULTS  TO  PROCESSES 


IN  AQUEOUS  SOLUTIONS 


5 . 1  Comparison  Between  the  Hydration  Energies  of  the 
Alkali  and  Halide  Ions. 

5.1a  In  Solution. 

The  heats  of  hydration  AHh  for  the  reaction  5.1  are 
obtained  from  the  heats  of  solution  at  infinite  dilution. 


5.1 


M+(g)  +  X  (g)  "  M+(ag)  +  X  (aq) 


o  o 

the  entropies  of  hydration  As  ^  from  S  and  the  Sackur- 

o  o 

Tetrode  equation  and  the  free  energies  AG  h  from  AH  h  and 

AS  ^ .  An  excellent  review  on  the  methods  of  assigning 

heats,  free  energies  and  entropies  to  single-ion  (M+^_^ 

or  X” ,  ,)  is  given  by  Rosseinsky  (31)  .  The  review  indi- 

(aq) 

cates  that  most  methods  are  based  on  either  dividing  by  2 
the  thermodynamic  functions  of  selected  anion-cation  pair 
or  using  some  extrapolation  procedure  involving  the  ionic 
radii.  Since  water  is  not  a  dimensionless  and  structur 
less  liquid,  the  hydration  of  equal  size  cations  and  anions 
is  not  necessarily  the  same  (  6 ) .  Latimer,  Pitzer  and 
Slansky  ( 32)  have  attempted  to  account  for  the  specificity 
in  ionic  hydration  by  adding  different  increments  to  the 
crystal  radii  of  anions  and  cations.  Verwey  (33)  has 
considered  the  different  orientations  of  water  molecules 


*• 
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around  ions  and  Buckingham  (34)  has  approached  the  problem 

by  taking  into  account  the  quadrupole  moment  of  the  water 

molecule.  The  agreement  between  these  different  methods 

is  often  surprisingly  good,  especially  in  the  case  of 

AH  ,  .  Sti31,  none  of  these  absolute  values  seem  entirely 
h 

reliable  since  they  are  all  based  on  the  assumption  that 
the  true  ionic  radii  are  known.  The  only  set  of  absolute 
hydration  energies  of  single  ions  M+  ^  ^  and  X  without 

any  assumption  about  the  dimensions  of  the  ions  is  the  one 
by  Randles  (35)  based  on  the  measurement  of  the  Volta 
potential  between  an  electrolyte  solution  and  another 
surface  separated  from  the  solution  by  vacuum  or  an  inert 

gas  . 

The  two  sets  of  hydration  energies  which  have  found 
wide  acceptance  are  those  of  Latimer,  Pitzer  and  Slansky 
and  Randles.  These  two  sets  were  used  in  the  previous 
discussion  3.2  for  a  correlation  of  the  relative  heats 
of  hydration  of  alkali  ions .  The  same  sets  of  data  will 
be  used  in  the  latter  part  of  this  discussion  for  the  cor¬ 
relation  of  the  present  results  with  the  total  single  ion 
hydration  energies. 


5  .  lb  In  the  Gas  Phase 

0 

The  gas  phase  hydration  enthalpies  AH  n_g/n  f°r  f^e 
alkali  (M+)  and  the  halide  (X  )  ions  are  shown  together 
in  Figure  5.1.  The  data  for  halide  ions  are  taken  from 
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positive  and  halide  negative  ions 
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the  studies  of  Arshadi  (30 )  which  were  made  in  this  lab¬ 
oratory.  The  hydration  interactions  of  these  spherical 
ions  are  probably  best  examined  if  one  compares  the  hyd- 
ration  energies  of  isoelectronic  pairs  like  K  and  Cl  , 

Rb+  and  Br~  etc.  Since  the  nuclear  charge  of  the  negative 
ion  in  the  isoelectronic  pair  is  two  units  lower,  the 
size  of  the  electron  cloud  of  the  isolated  negative  ion 
will  be  definitely  larger  than  that  of  the  positive  ion. 
Therefore,  if  one  considers  only  the  effect  of  ionic  oize 
(in  vacuum)  one  would  definitely  predict  that  the  posiuive 
isoelectronic  ion  should  have  higher  -]-n  Chapter 

4 . 7  we  have  also  suggested  that  the  initial  steps  of  hyd¬ 
ration  of  positive  ions  are  assisted  by  the  formation  of 
a  dative  bond  between  the  oxygen  atom  of  the  water  molecule 
and  the  ion.  The  existence  of  such  a  dative  bond  between 
isoelectronic  negative  ions  and  water  molecules  is  not 
possible.  As  explained  in  a  previous  Chapter  (4.7)  the 
effect  of  the  dative  bond  decreases  with  increasing  size 
of  the  positive  ion.  If  one  examines  the  enthalpies  of 
the  larger  ions  like  Rb+  and  Br“,  Cs+  and  l“,  for  which 
the  effect  of  any  covalent  dative  bonding  to  water  should 
be  small,  one  still  notices  (Figure  5.1)  that  isoelect- 

O 

ronic  positive  ions  have  considerably  higher  “AHn„1/n* 

Thus  the  present  results  of  hydration  energies  of  iso¬ 
electronic  positive  and  negative  ions,  support  at  least 
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qualitatively  the  prediction  based  only  on  the  size  of 
the  ions.  It  is  interesting  to  note  that  the  difference 

O 

between  the  AH  ,  of  the  isoelectronic  pairs  decreases 

n-l,n 

+  _  - 

as  n  becomes  larger  and  that  in  the  cases  of  Na  and  F  , 

K+  and  Cl  ,  Rb+  and  Br  a  crossover  is  observed  within 

O 

the  range  of  the  figure  in  which  the  -AH  ,  of  the 

1 1  -L  /II 

negative  ion  becomes  higher.  Furthermore  the  smaller  the 
size  of  the  ions,  the  sooner  the  crossover  occurs.  For 
Na"1  ,  F  crossover  occurs  at  about  n  =  4,  for  K+ ,  Cl  at 
n  *  5.5  and  for  Rb+,  Br”  at  n  :  6.  These  observations 
can  be  explained  on  the  basis  of  electrostatic  calculations 
which  indicated  (30)  that  the  slower  decrease  of  “AHn_1^n 
of  the  negative  ions  is  due  to  the  ability  of  the  water 
molecules  to  assume  nonsymmetric  positions  in  which  only 
one  of  the  hydrogen  atoms  "touches"  the  negative  ion.  It 
has  been  shown  (30)  that  these  nonsymmetric  arrangements 
result  in  lower  water-water  repulsions  for  ligands  in 
the  same  shell  and  become  of  importance  when  the  shell 
becomes  crowded.  Of  course,  crowding  occurs  soon  in  the 
smaller  ions. 

It  is  of  interest  to  examine  at  this  point  the  cor¬ 
relation  of  the  present  results  with  the  total  single  ion 
hydration  enthalpies.  Values  for  the  total  enthalpies 
of  hydration  have  been  obtained  by  several  authors  dis¬ 
cussed  previously  (31) .  We  used  two  sets  of  data  for 
total  enthalpies  of  hydration,  those  of  Latimer,  Pitzer 
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and  Slansky  (32)  and  Randles  (35).  The  same  sets  were 
used  earlier  in  Chapter  (  3  )  for  correlation  of  the 
relative  heats  of  hydration  of  alkali  ions . 

Examining  the  values  that  the  Latimer  set  gives 
for  the  alkali  and  halide  ions  one  finds  that  the  -AH^ 
for  the  negative  ion  is  substantially  larger  than  that 
for  the  positive  ions.  Thus  -AH^(Br  )  =  81.4  kcal/mole 
while  -AHh(Rb+)  =  69.2  kcal/mole.  Therefore,  Latimer's 
heats  of  hydration  are  exactly  of  opposite  magnitudes 
to  those  expected  from  the  standpoint  of  electron  cloud 
size  of  the  gaseous  alkali  and  halide  ions.  Two  basi¬ 
cally  different  reasons  have  been  put  forward  to  explain 
the  Latimer  data.  The  first  explanation  which  exists  in 
several  variations  (32 )  proposes  that  the  arrangement  of 
water  molecules  around  the  negative  ion  and/or  in  the 
transition  from  the  hydrated  ion  to  the  bulk  liquid  is 
more  favorable.  The  second  suggestion  has  been  that  water 
molecules  have  an  electrical  quadrupole  moment  which  is 
of  such  a  size  that  it  leads  to  an  attraction  with  neg¬ 
ative  ions  and  repulsion  with  positive  ions  (34) -  On 
the  basis  of  the  present  results  which  give  higher 
energies  of  interactions  for  the  positive  ions,  it  is 
possible  to  eliminate  the  quadrupole  theory  since  the 
0ff0ct  of  the  quadrupole  should  be  particularly  large  at 
close  range  i.e.  exactly  the  condition  present  in  the  gaseous 
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clusters.  It  should  be  pointed  out  that  the  present  re¬ 
sults  do  not  prove  that  a  quadrupole  moment  is  completely 
absent.  They  only  show  that  a  water  quadrupole  moment 
of  the  magnitude  required  to  explain  the  Latimer  data 
cannot  be  present.  In  any  event,  as  is  shown  in  the 
subsequent  discussion,  the  present  results  cannot  be 
reconciled  with  the  very  strong  hydration  interactions 
for  negative  ions  predicted  by  the  Latimer  data. 

The  Randles  total  hydration  energy  set  assigns  a 
slightly  higher  -AHh  for  F~  than  for  Na+  but  for  the 
larger  ions  the  — AH^  of  the  positive  isoelectronic  ions 
are  somewhat  larger  than  those  of  the  negative  ions .  The 


Randles  data  are  thus  closer  to  expectations  based  on 
the  present  results .  A  correlation  of  the  present  data 
and  the  single  ion  hydration  energies  of  Latimer  and 


Randles  is  given  in  the  composite  Figure  5.2a,b  which 
shows  plots  of  [-AH  (M  ) ]  -  [-AH  (X  ) ]  for  isoelect- 

ronic  pairs.  The  plots  are  based  on  the  experimental 

o 

AH  whenever  these  were  available  (low  n) .  For 

n-l,n 

higher  n  the  values  from  Figure  5.1  were  used.  Also 
indicated  in  Figure  5.2  (as  two  horizontal  lines)  are  the 
differences  of  the  total  single  ion  hydration  energies 
for  the  isoelectronic  pairs  based  on  the  data  of  Latimer 
and  Randles.  It  can  be  seen  that  the  AHQ^n  differences 
are  in  all  cases  much  closer  to  the  Randles  differences, 
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FIGURE  5a  Comparison  of  AH  hydration  of  alkali  and  halide  ions 
with  total  enthalpies  of  hydration  of  Latimer,  Pitzer 
and  Slansky  (L)  and  Randles  (R) . 
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Slansky  (L)  and  Randles  (R) . 
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moreover  the  AH  differences  seem  to  extrapolate  for  high 

vJ  j  1 1 

n  into  the  Randles  differences.  Of  course  by  definition 
the  AHq  ^  differences  should  become  equal  to  the  AH^ 
differences  for  very  high  n. 

We  can  summarize  the  observations  on  positive  and 
negative  isoelectronic  pairs  as  follows .  In  the  initial 
hydration  interactions  (n  small)  the  positive  ion  gives 
higher  energies  of  interaction  as  would  be  expected  from 
its  smaller  size.  As  the  cluster  grows  interactions  in 
the  negative  ion  become  gradually  more  favourable  (Figure 
5.1  and  5.2) .  This  is  probably  due  to  the  ability  of  the 
negative  ion  to  pack  the  water  molecules  more  closely 
without  too  large  water-water  repulsions.  Extrapolation 
of  the  data  to  moderately  high  n  seems  to  lead  to  an 
asymptotic  approach  to  the  Randles  data.  Thus  the  present 
results  are  in  support  of  the  Randles  data  and  also  indi¬ 
cate  that  the  decisive  interactions  determining  the  total 
heats  of  hydration  occur  during  the  attachment  of  the  first 
8-12  molecules.  One  may  consider  the  further  path  towards 
formation  of  a  liquid  solution  as  a  single  step  in  which 
the  large  clusters  are  fitted  into  the  bulk  of  the  liquid. 
The  present  results  indicate  that  the  hydration  energy  for 
this  final  "step"  will  be  quite  similar  for  a  positive  or 
negative  cluster  containing  the  corresponding  isoeletronic 


ion . 
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5b.  Ion-Ligand  Bond  Energies  and  Hydration  Numbers  in 

Solution . 

According  to  Bockris  (  5  )  the  term  hydration  number 
when  applied  to  simple  ions  in  solution  corresponds  to  the 
coordination  number  of  the  ion  in  solution.  Since  the 
water  molecules  coordinated  to  the  ion  are  in  rapid  exchange 
with  the  bulk  of  the  solution  it  is  only  on  a  time  average 
that  we  may  speak  of  some  molecules  being  immobilized  by 
the  ionic  field.  The  hydration  number  is  therefore  an 
empirical  parameter  which  gives  the  effective  number  of 
water  molecules  that  have  undergone  some  constant  critical 
change  in  property  related  to  a  suitable  definition  of  hyd¬ 
ration  (  6 ) .  It  has  been  suggested  (  6 )  that  whenever 
possible,  properties  at  infinite  dilution  should  be  used 
for  the  determination  of  hydration  numbers  in  order  to 
avoid  any  overlap  with  other  types  of  interactions  in  the 
solutions.  Since  the  hydration  in  the  gas  phase  may  be 
compared  to  the  properties  of  hydration  at  "infinite 
dilution"  in  solution,  one  way  of  evaluating  the  hydration 
numbers  of  alkali  ions  would  be  the  use  of  present  data. 

O 

One  way  in  which  the  present  AH  n_1  n  values  can  be 
related  to  the  liquid  phase  is  through  the  theoretical 
arguments  of  Bernal  and  Fowler  (77) •  They  have  proposed 
the  following  relating  to  the  ionic  solvations  in  liquid 
water.  The  hydration  of  an  ion  will  occur  when  the 
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potential  energy  of  a  water  molecule  coordinated  to  an  ion 
is  less  than  that  of  a  free  water  molecule.  The  free  water 
molecule  has  four  neighbors  while  a  single  water  molecule 
bonded  to  an  ion  has  three/  namely  the  ion  and  two  water 
molecules.  The  condition  for  an  ion  having  at  least  one 
coordinated  water  molecule  is,  therefore,  that  the  potential 
energy  AE  n  ^  n(l)  ^or  a  water  ^dc011!6  due  to  coordination 
to  an  ion  be  less  than  that  due  to  two  other  water  mole¬ 
cules.  The  electrostatic  potential  energy  due  to  two 
water  molecules  is  calculated  by  Bernal  and  Fowler  to  be 
-15  kcal/mole.  This  fact  gives  rise  to  equation  5.2  which 
expresses  the  ion— water  molecule  bond  energy  B.L. 


5.2 


B-E-  =  &E  n-l , n ( 1)  +  15 


AE°n-l,n(l)  can  be  SSt  eqUal  t0  AE  n-l,n(g)'  the  Potent:La:L 
energy  for  a  water  molecule  due  to  coordination  to  an  ion 

in  the  gas  phase.  This  is  not  unreasonable  since  Bernal  and 

Fowler  appear  to  have  taken  the  dielectric  constant  as  unity 

in  their  calculation  of  the  electrostatic  potential  energy 

due  to  two  water  molecules. 

o 

Since  in  section  4.7  it  was  found  that  AH  n_g^n 

°  R°  ~  AE°  ,  _m  Equation  5.2  can  be 

AE  n-l ,n  (g)  n-l,n  n-l,n(l) 

rewritten  as  equation  5 . 3 


5.3 


r  F  =  AH  ,  +15 

B.-h.  an  n_x^n 
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The  bond  energies  calculated  for  M  (H^0^  -  H^O  with  the 

Ah  values  from  Tables  3.1  to  3.4  are  listed  in 

n-1 ,  n 

Table  5.1.  The  bond  energies  for  M+(H20)n  -  H20  show 
that  the  Li+  in  solution  may  have  up  to  4  coordinated 
water  molecules  (hydration  number  ~4) ,  Na  ion  up  to  3, 

I  _L,  -{- 

K  ion  up  to  2 ,  Rb  ion  may  have  one  and  Cs  ion  no  coor¬ 
dinated  water  molecules,  (hydration  number  ~0) .  The 
hydration  numbers  obtained  show  a  regular  decrease  with 
increasing  size  of  the  ion.  For  Li  to  Cs  the  hydration 
number  varies  from  4  to  0.  It  has  been  pointed  out  (6) 
that  ions  of  low  charge  and  large  radius  (Cs  for  example) 
cannot  form  ion— ligand  bonds  in  solution  (hydration 
number  ~0)  ,  since  there  will  be  a  critical  radius  above 

which  the  charge  density  falls  below  the  value  required 
to  reorient  water  molecules.  At  larger  values  of  the 
radius  of  the  ion  there  is  no  preferred  orientation  of  the 
water  molecules.  Thermal  agitation  is  then  able  to  breah  up 
the  water  structure  in  the  vicinity  of  the  ion.  The  original 
calculations  of  Bernal  and  Fowler  (77)  indicated  that  for 
univalent  ions  this  critical  radius  would  be  about  1.6A. 
According  to  Pauling  (69) ,  Goldschmidt  (61)  and  Gourary 
and  Adrian  (70)  the  ionic  radii  of  Cs  are  1.69,  1.65  and 
1.83  respectively.  Thus  if  the  original  calculations  of 
Bernal  and  Fowler  are  correct  one  would  expect  a  zero 
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TABLE  5 . 1 

I  on -water  Molecule  Bond.  Energies  in  the  Liquid. 


Phase 

Bond 

Bond  energy 
(kcal/mole) 

Li+  H20 

-19.0 

Li+  (H20)  — 

H2° 

00 

• 

o 

1 — 1 

1 

Li+  (H20) 2 — - 

'  H2° 

-  5.7 

Li+  (H20)  3 — - 

'  H2° 

-  1.4 

Li+(H20)  4— 

'  H2° 

+  1.1 

Na+ - H20 

-  9.0 

Na+(I420)  - 

-  H2° 

-  4.8 

Na+  (H20) 2  — 

H2° 

-0.8 

Na  (H20)  3 

h2o 

+  1.2 

K+  H20 

-  2.9 

k+(h20)  - - 

o 

CM 

iU 

-  1.1 

K+(H20)2 

-  h20 

+  1.8 

Rb*  - —  H20 

-  0.9 

Rb+(H20)  - 

h20 

+  1.4 

Cs+  H20 

+  1.3 

144. 


TABLE  5.2 


Hydration  Numbers  h  for  Alkali  Ions 


Ion 

a 

Hydration  numbers  h 

b  c 

d 

Li+ 

4 

4.2-4. 7 

4 

4 

Na+ 

3 

2 . 7-3 . 3 

3 

3.1 

K+ 

2 

CO 

• 

1 - 1 

1 

r- 

* 

i — l 

2 

2.1 

Rb+ 

1 

1.5 

1 

1.6 

Cs  + 

0 

• 

1 - 1 

0 

o 

• 

i — 1 

(a)  Present  calculation 

(b)  Modified  activity  coefficient  equations, 

ref .  (78 )  . 

(c)  Mean  values  from  methods,  mobility,  entropy, 
compressibility,  density,  ref.  (  5  ) • 

N.M.R.  effects  in  aqueous  solutions  of  electro¬ 
lytes  ,  ref .  (79  )  . 


(d) 
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+ 

hydration  number  for  the  Cs  ion. 

The  values  of  hydration  numbers  calculated  here  are 
very  close  to  those  obtained  from  other  types  of  measure¬ 
ments,  described  in  references  (78)  (5  )  (79).  See  Table 

5.2. 
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6 .  OTHER  APPLICATIONS  OF  THE  GAS  PHASE  HYDRATION 

MEASUREMENTS 

6 . 1  Clustering  Equilibria  and  the  Catalysis  of  Water 

Condensation  by  Ionic  Charges. 

O 

It  can  be  seen  from  Figure  3.27  where  “AGn_^  n  is 
plotted  versus  n,  that  with  increasing  n  (n  >>  6)  the 

O 

-AG  _  values  will  approach  the  free  energy  of  evapora- 
n-l,n 

tion  -AG  /ooootvm  =  2.03  kcal/mole  for  reaction 
evap (29  8  °K) 


H2°  (g)  H2°(St) 

In  other  words  as  the  ionic  charge  in  the  cluster  M  (^2^  n 
becomes  more  and  more  shielded  by  water  molecules  (large 
increase  of  n)  the  cluster  M  (H^O)  will  approach  a  liquid 
droplet . 

In  Figure  6.1  the  -AG°  ,  values  at  298°K  were 

II  ~L  f  1 1 

plotted  versus  n,  using  a  logarithmic  scale  for  n  in  order 
to  be  able  to  accomodate  a  large  increase  of  n.  The  dashed 
line  is  calculated  by  equation  6.1 


6.1 


O 


AG 


evap 


o 


AG 


(2) 


where  the  free  energy  change  AG^j  is  equivalent  to  the 
transfer  (of  one  mole  of  water)  from  pressure  P  —  23.7  torr 
(water  vapour  pressure  at  298°K)  to  P  =  1  atm,  given  by 
equation  6 . 2 


-RT 


In 


760 


6.2 


23.7 


FIGURE  6.1 

O 

Comparison  of  free  energy  of  hydration  hGn_1^n 
with  free  energy  of  evaporation  of  1  mole  of 
water  molecules  from  1  mole  of  water  droplets 
carrying  a  single  ionic  charge.  Number  of 
molecules  contained  per  droplet  given  in  the 
abscissa.  Free  energy  of  evaporation  calculated 
from  J.  J.  Thomson  formula. 
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FIGURE  6.1 


Number  of  Water  Molecules 


148. 


The  AG^-jj  represents  normal  free  energy  of  evaporation 
from  a  flat  water  surface  at  298°K. 

The  AG^)  gives  the  free  energy  of  evaporation  of 
water  from  the  spherical  droplet  of  radius  r,  carrying  a 
single  centrally  located  charge.  The  AG^j  values  as  a 
function  of  radius  of  the  droplet  are  represented  by 
the  J.  J.  Thomson  formula  (80)  equation  6.3 


6.3 


AG 


(2) 


RT  In 


P  r 


8TTer 


Pr  is  the  vapour  pressure  over  the  droplet  while  P^is  the 
vapour  pressure  of  water  with  a  flat  surface.  p  and  M  are 
the  density  and  molecular  weight,  while  0,  e  and  e  are  the 
surface  tension,  electronic  charge  and  dielectric  constant 
The  equation  6.3  gives  the  condition,  when  a  droplet  of 
radius  r  will  be  in  equilibrium  with  the  surrounding 
medium,  that  is  the  droplet  will  neither  shrink  nor  grow. 
The  first  term  (~)  in  the  equation  represents  the  lower- 

O 

ing  of  the  AG  with  decreasing  radius  of  the  droplet. 

This  effect  is  due  to  the  operation  of  surface  tension. 

The  vapour  pressure  of  the  spherical  droplet  is  greater 
than  that  of  the  flat  surface,  the  pressure  increasing  as 


the  radius  of  the  droplets  decreases.  The  higher  vapour 
pressure  of  small  droplets  accounts  for  their  tendency 
to  disappear  by  "distillation"  to  larger  droplets.  Large 
droplets  thus  tend  to  grow  at  the  expense  of  the  smaller 
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2 

0 

ones.  The  second  term  (  ■  -  -  -^-)  in  the  equation  6.3  expresses 

8Tier 

the  increase  of  AG^)  (becomes  less  negative)  with  the 
decreasing  radius  of  the  droplet.  This  effect  is  due  to 
the  increased  attraction  between  the  charge  and  the  water 
dipoles .  As  the  radius  decreases  the  distinction  between 
charged  water  droplets  and  water  clusters  m"1  (^O)  becomes 
less  pronounced.  The  P^/P^  and  consequently  the  AG^) 
approach  zero  values ,  due  to  the  increased  attraction 
between  the  charge  and  water  dipoles  (second  term  in 
equation  6.3}.  Thus  the  AG  evap  by  equation  6.1  begins  to 

o 

approach  AG  . 

n  _L  f  li 

In  Figure  6.2,  taken  from  reference  (80),  the 
variation  of  ^r/P  as  a  function  of  droplet  radius  is 
plotted  as  a  broken  line  for  the  case  of  water  vapour  and 
as  a  solid  line  for  a  singly  charged  water  vapour  droplet 

O  o 


at  29  8  °K .  Tab-le  6.1  gives  the  values  of  AG  ^  ,  AG  ^ 

and  AG  as  a  function  of  the  radius  of  the  charged 

evap 

droplet.  The  AG  values  were  obtained  with  macroscopic 

values  for  p ,  0  and  e . 


p 

The  calculated  values  for  -  as  a  function  of  the 

00 

radius  of  the  droplet,  given  in  Table  6.1,  are  compared 

_  _  g 

with  the  values  from  Figure  6.2.  For  the  r  >  7  x  10  cm, 

P 

the  calculated  r /Pm  values  are  consistent  with  those  from 
Figure  6.2.  However  for  r  <  7  x  10  ,  the  r/P00  ratios  were 

taken  from  Figure  6.2  (solid  line),  since  the  Pr/P0O  values 
for  very  small  droplets  cannot  be  calculated  with  macro- 
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scopic  values  for  p,  0  and  e  (81) . 

Figure  6.1  shows  that  the  experimentally  determined 

O 

AG  ,  values  for  all  ions  are  consistent  with  the  pre- 
n-l,n  c 

O 

dictions  of  equation  6.3  since  the  AG  n  for  all  ions 

extrapolate  smoothly  for  high  n  into  the  calculated  curve. 
The  Aq  ,  values  for  H'  (Ho0)  clusters,  presented  in 
this  figure  were  taken  from  reference  (24). 

Figure  6.1  also  illustrates  the  fact  that  there  is 
a  theoretical  limit  for  n  beyond  which  clustering  equilib- 

O 

ria  cannot  be  measured.  Thus  when  Ag  becomes  smaller 

n  ~L  y  II 

than  AG  =  2.063  kcal/mole  at  298°K  (normal  free  energy 

evap 

of  evaporation  from  flat  water  surface  indicated  by  a  black 
horizontal  line  in  Figure  6.1)  the  condensation  to  bulk 
liquid  water  becomes  theoretically  possible.  However,  in 
the  absence  of  other  nucleation  catalysts,  condensation 
will  not  occur  and  ion  clustering  equilibria  will  exist  (at 
water  pressure  above  saturation  pressure)  until  n  becomes 

O 

about  equal  to  20.  Beyond  this  point  the  Ag  n_1  n  (dashed 
curve)  begins  to  increase  which  means  that  clusters  with 
higher  n  become  thermodynamically  more  stable  and  growth 
of  the  clusters  becomes  spontaneous  ending  in  bulk  water. 

In  experiments  done  up  to  the  present  we  found  that 
clustering  equilibria  could  not  be  measured  beyond  n  ~  7 
due  to  condensation  of  water  in  the  ion  source.  According 
to  the  above  discussion  the  condensation  must  have  been 


. 
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catalyzed  by  some  other  nucleating  agents  (i.e.  active 
points  on  the  metal  surfaces  of  the  walls  etc.).  Thus 
application  of  special  techniques,  for  instance  those 
used  in  cloud  chamber  experiments  coupled  with  rapid  mass 
scanning  might  allow  equilibria  to  be  measured  for  n  up 
to  about  20,  but  not  higher. 

6.2  Proton  Affinities  (Basicities)  of  the  Alkali  Hyd¬ 
roxides  in  the  Gas  Phase. 

One  interesting  result  from  the  present  work  is 
that  enthalpy  changes  Ahq  for  the  reaction  M+  +  H^O 
(MOH2)+  (0,1)  when  combined  with  available  thermochemical 

data  allow  the  determination  of  the  proton  affinities  of 
the  (gaseous)  alkali  hydroxides.  These  proton  affinities 
are  of  considerable  chemical  interest  especially  when 
compared  to  those  of  other  related  compounds. 

The  proton  affinity  (P.A.)  of  the  hydroxides,  MOH 
is  defined  as  the  negative  of  the  enthalpy  change  for  the 
gas  phase  reaction  6 . 1 

6.1  MOH  +  H+  (MOH2)+  P.A.  (MOH)  =  -AH 

The  enthalpy  change  for  reaction  6.1  can  be  evaluated  from 
the  enthalpy  of  formation  of  the  reactants  and  the  product. 
The  enthalpies  of  formation  of  MOH  and  H+  are  available. 
The  enthalpy  of  formation  of  (MOH^4  can  be  evaluated 
from  reaction  (0,1)  since  AHQ  1  was  determined  and  AHf(M+) 
and  AH  (HO)  are  available.  Thus  the  equation  6.4 
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6.4  AHq  x  =  AHf(MOH2)+  -  AHfH20  -  AHfM+ 

when  combined  with  the  equation  6.5  representing  the  en¬ 
thalpy  change  for  reaction  6.1. 

6.5  AH  =  AHf(MOH2)+  -  AHfMOH  -  AHfH+ 

would  give  the  proton  affinity  of  MOH  represented  by 
equation  6 . 6 

6.6  P . A . (MOH)  =  AHf (MOH)  -  AHQ  ±  -  AHf(M+)  +  AHf(H+) 

AHf  (H20) 

The  results  and  the  data  used  for  their  calculation  are 
given  in  Table  6.2.  The  values  AHf(Ii+)  =  367  and  AHf(H20) 
-57.8  kcal/mole  were  also  used.  The  proton  affinities  of 
the  hydroxides  which  are  in  the  range  240  -  270  kcal/mole 
are  considerably  larger  than  those  of  water,  methanol  and 
ammonia  which  are  approximately  160,  185  and  206  kcal/mole 
respectively  (82) .  It  has  been  pointed  out  (83)  that  in 
the  series  R— OH  the  proton  affinity  should  increase  with 
the  electron  donating  ability  of  the  substituent  R  (84). 
The  very  high  proton  affinity  of  the  alkali  hydroxides  can 
be  thus  understood  in  terms  of  the  very  high  electron 
donating  ability  of  the  substituent  M.  This  ability  is 
so  extreme  that  in  the  protonated  species  essentially 
all  the  charge  resides  in  M  i.e.  the  resulting  structure 
is  that  of  the  hydrated  alkali  ion  as  indicated  by  the 
conventional  formula  M  OH2 . 
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The  variation  of  the  proton  affinity  in  the  MOH  series 
with  the  electron  releasing  ability  of  M  is  illustrated  in 
Figure  6.3f  where  the  proton  affinities  are  plotted  as  a 
function  of  the  ionization  potential  of  the  alkali  atoms. 

The  rationale  of  the  plot  being  that  the  electron  releasing 
ability  decreases  as  the  ionization  potential  of  the  alkali 
ion  increases.  It  can  be  seen  that  an  almost  linear  relation¬ 
ship  is  obtained.  Considering  the  proton  affinity  as  a 
measure  of  the  base  strength  of  a  (gas  phase)  Bronsted 
base  we  find  that  the  base  strength  of  hydroxides  increases 
from  LiOH  to  CsOH.  This  order  of  base  strength  need  not 
be  observed  in  solutions  since  there  the  process  is  com¬ 
plicated  by  solvation  effects. 

Since  the  proton  affinity  (P.A.)  of  a  hydroxide  is 
defined  as  the  negative  of  the  enthalpy  change  for  the 
reaction  6.1f  the  hydrogen  affinity  (II. A.)  can  be  simil¬ 
arly  defined  as  the  negative  of  the  enthalpy  change  for 
the  reaction  6.2. 

6.2  MOH+  +  H  -*  (MOH2)  + 

6.7  H .  A .  (MOH+)  =  AHf(MOIl+)  +  AHf  (Ii)  -  AHf(MOH2)  + 

Relating  P.A.  (MOH)  to  H.A.(M0I1+),  the  following  expression 
represented  by  equation  6.8  is  obtained 


6.8  P.A. (MOH)  -  H.A.(MOH+)  =  I.P.(H)  -  I.P.(MOH) 
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where  I.P.(H)  and  I.P.(MOH)  are  the  ionization  potentials 
of  H  and  MOH  respectively.  Thus  the  proton  and  hydrogen 
affinities  of  hydroxides  and  their  corresponding  ions  are 
interrelated  by  ionization  potentials  as  shown  in  the 
above  equation  6.8.  The  hydrogen  affinity  H.A.(MOH+)  is 
simply  the  MOH+-  H  bond  strength  D (MOH+-  H) .  From  equation 
6.7  D(NaOH+-  H)  =  140  kcal/mole  was  calculated.  The  values 
I.P.(H)  =  313  and  I.P.(NaOH)  =  207  kcal/mole  (85)  were 
used.  Unfortunately  the  ionization  potentials  for  other 
alkali  hydroxides  are  not  available  in  the  literature 
therefore  the’ hydrogen  affinities  for  other  alkali  hyd¬ 
roxides  were  not  calculated. 

The  hydrogen  affinity  H.A.  (NaOH+)  =  140  kcal/mole, 
can  be  compared  to  that  of  water,  H.A.  (H20+)  =  141 
kcal/mole  (86) .  Assuming  that  the  values  of  P.A.  (NaOH  = 

247.6  and  P.A.(H20)  =  164  kcal/mole  are  of  right  order  of 
magnitude,  the  difference  in  proton  affinities  AP.A.  = 

83.6  kcal/mole  would  be  compensated  by  AI.P.  =  83 
kcal/mole,  since  I.P.(H20)  =  290.6  and  I. P.  (NaOH)  =  207 
kcal/mole.  Thus  according  to  equation  6.7  D (HOH  -H)  ~ 
D(NaOH+-H)  ~  140  kcal/mole. 
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7  *  SUGGESTION  FOR  FURTHER  RESEARCH 

7.1  '  Solvation  of  Alkali  and  Halide  Ions  by  Polar  Aprotic 

Solvents . 

Since  the  gas  phase  solvation  of  alkali  ions  plays  an 
important  role  in  gas  phase  chemistry  generally  and  has 
particular  relevance  to  the  solvation  in  solutions,  it  is 
proposed  that  further  work  be  carried  out  in  this  area. 

The  energetics  of  ion  solvation  in  nonaqueous  solutions 
have  recently  become  the  subject  of  many  investigations, 
for  example  (87,88).  Considerable  evidence  indicates  (89) 
that  in  polar  aprotic  solvents,  such  as  acetonitrile 
(CH^CN) ,  acetone  (CH^COCH^) ,  nitromethane  (CH2N02)  etc., 
negative  ions  are  solvated  more  weakly  than  positive  ions. 

It  was  also  observed  that  all  negative  ions  are  solvated 
more  weakly  in  aprotic  solvents  than  in  water  (89) . 

Although  the  dipole  moments  of  aprotic  solvents  are 
much  higher,  acetonitrile  for  example  has  a  dipole  moment 
of  3.37D,  compared  to  1.85D  for  water,  water  is  still  a 
better  solvating  agent  for  negative  ions  than  acetonitrile. 

It  has  been  suggested  that  this  effect  is  due  to  steric 
inhibition  (89) .  Because  of  the  vicinity  of  the  electron 
cloud  of  the  C=N  bond,  an  anion  probably  will  be  forced 
to  take  up  the  specific  position  indicated  in  structure  I, 
resulting  in  loss  of  freedom  of  rotation.  Furthermore, 
solvation  of  a  (small)  anion  by  more  than  one  acetonitrile 
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molecule  then  becomes,  improbable.  The  same  is  not  true  for 
the  cylindrically  symmetrical  solvation  of  a  positive  ion 
indicated  in  structure  XI. 


H 


H - C - CEN 


/ 

X 


E 

H - C - CEN - M+ 

I 

E 


I 


II 


Therefore  it  would  be  very  desirable  to  measure  the  gas 
phase  solvation  of  positive  alkali  and  negative  halide 
ions  by  acetonitrile.  The  data  obtained  would  indicate 
how  many  molecules  of  acetonitrile  are  attached  to  each 
particular  ion  in  the  first  and  second  shell.  Furthermore , 
the  enthalpies  for  the  positive  and  negative  ions  will 
show  directly  whether  the  solvation  by  negative  ions  is 


much  weaker. 
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